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CHAPTER  1 
INTRODUCTION 

Gyrotrons  belong  to  a  class  of  vacuum  electron  devices 
whose  operating  principle  is  based  on  the  relativistic 
interactions  between  the  transverse  rf  fields  and  the  gyrat¬ 
ing  electron  beams,  [l]-[8].  The  basic  physical  principle 
underlying  the  operation  of  gyrotron  devices  is  sometimes 
referred  to  as  cyclotron  maser  instability  (CMI).  It  is 
generally  aknowledged  that  Twiss  [9]  in  1958  was  the  first 
to  recognize  it  as  an  amplifying  mechanism  for  free  electron 
gyroradiat ion .  Since  then,  research  in  the  CMI/gyrotron 
effect  has  been  very  active  in  understanding  the  phenomenon, 
finding  applications  and  tapping  this  new  source  of  radia¬ 
tion  . 

Most  of  the  work  on  gyrotrons  has  concentrated  in  the 
past  in  achieving  high  power,  efficient  MM  wave  oscillators 
for  plasma  heating  experiments.  Considerable  success  has 
been  reported  in  this  area.  Andronov  et  al,[10],  reported  on 
experiments  with  gyrotrons  that  operated  at  MM  wavelength 
with  output  power  over  1MW  within  a  pulse  duration  of  about 
10~*s.  In  one  of  the  experimental  gyrotrons  operating  at  3mm 
wavelength,  a  pulsed  ouput  power  of  as  much  as  2  .  1MW  was 
reportedly  attained,  [11]  cited  in  [15].  Other  frequencies 
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were  investigated  and,  in  an  experiment  with  a  gyromonotron 
using  a  TE01  circular  mode  and  operating  at  350Hz,  a  peak 
output  power  of  about  150kW  with  a  31%  efficiency  at  lOOkW 
was  achieved,  [12];  a  gyrotron  operating  at  the  same  35GHz 
frequency  but  using  a  TEQ4  mode,  [13],  achieved  output 
powers  of  320kW  with  an  efficiency  of  40%.  Boucher  et  al  . 
[14],  achieved  150kW  output  power  at  the  same  frequency 
( 3  5GHz )  using  a  TE 02  circular  mode.  Fix,  Flyagin  et  al  . 
reported  in  [15]  results  such  as  lOOkW  quasi~CW  output  power 
of  3mm  wavelength,  150kW  with  10%  efficiency  at  5.6mm,  140kW 
with  25%  efficiency  at  9mm.  Jory,  Friedlander  et  al . ,  [16], 
built  a  pulsed  gyrotron  which  produced  248kW  peak  output 
power  at  28GHz  and  34%  efficiency,  Symons  and  Jory,  [17], 
achieved  212kW  CW  output  power  at  28GHz. 

Although  gyrotrons  oscillators  have  had  some  success  as 
rf  sources  for  plasma  heating  in  fusion  devices,  their  lack 
of  signal  coherence  and  the  higher  signal  quality  required 
for  applications  in  radar  and  communications  has  eliminated 
them  as  potential  sources  in  those  applications.  However, 
another  possible  application  of  the  CMI  effect  has  been 
identified:  since  as  an  oscillator  it  is  not  a  source 
suitable  for  radar/communication  purposes,  it  is  of  interest 
to  look  into  the  possibility  of  achieving  high  quality 
amplifier  application  of  the  basic  CMI  effect. 
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To  be  regarded  as  a  high  quality  amplifier,  it  is  well 
accepted  that  the  device  should  have  large  instantaneous 
bandwidth,  and  also  nearly  linear  phase  versus  frequency 
characteristics.  It  is  in  these  aspects  that  the  gyrotron 
amplifier  effort  has  been  less  than  successful.  A  C-band 
gyro-TWT  experiment  at  Varian  was  performed  by  Symons,  Jory, 
Ferguson  and  Valier,  [18]  —  [19]  .  In  this  experiment,  the 
fundamental  electron  cyclotron  frequency  of  the  beam  inter¬ 
acts  with  the  TE^j  dominant  waveguide  mode:  an  output  power 
of  50kW  with  16.6%  efficiency  and  6%  bandwidth  was  achieved, 
[15].  Higher  power  outputs  or  larger  bandwidths  were  achie¬ 
ved  using  various  beam  powers.  The  largest  instantaneous 
bandwidth  was  9.3%  for  18.8kW  peak  output  power  at  a  9.8% 
efficiency  for  a  40kV-4A  beam  power.  Another  experiment  by 
Park,  Chu  et  al . ,  [20],  reported  a  20kW  power  output  at  6GHz 
with  a  10%  efficiency,  but  a  bandwidth  no  better  than  3.7%. 
Very  little  information  on  the  phase  characteristic  was 
reported  in  these  experiments. 

The  reason  for  the  narrow  bandwidth  is  well  understood: 
the  rapid  increase  of  the  group  velocity  as  a  function  of 
frequency  in  a  smooth  walled  waveguide  circuit  allows  the 
interaction  between  the  beam  and  the  wave  to  take  place  over 
only  a  very  narrow  band.  Chapter  2  shows  that  the  best 
operating  condition  is  when  the  beam  mode  and  waveguide  mode 
are  tangent,  the  so-called  grazing  condition.  In  this  case, 
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the  frequency  band  where  the  instability  is  excited  (i.e., 
where  amplification  can  occur)  is  centered  at  the  tangential 
point,  with  not  much  "elbow  room”  (see  figure  1.1).  Research 
to  overcome  this  drawback  is  very  active,  and  several  sche- 
mes  to  improve  the  bandwidth  have  been  suggested.  Barnett, 
Lau ,  Chu  and  Granatstein,  [21],  used  a  tapered  waveguide 
circuit  in  an  experiment  that  achieved  a  small  signal  3-dB 
bandwidth  of  13*  with  18-dB  gain  at  midband(35  GHz).  The 
same  investigators  suggested  another  structure  using  a 
tapered  waveguide  as  well  as  a  tapered  external  magnetic 
field  profile:  they  calculated  that  a  saturated  bandwidth 
greater  than  two  octaves  is  possible  for  both  TE01  and  TE ^  ^ 
modes,  [22].  A  tapered  two-stage  structure  with  a  tapered 
magnetic  field  was  proposed  by  Ganguly  and  Ahn,  [23]:  their 
calculations  showed  that  a  45*  bandwidth  at  45dB  gain  and 
25*  efficiency  is  possible  if  the  axial  velocity  spread  of 
the  beam  is  zero;  they  found  however  that  both  the  gain  and 
the  bandwidth  are  very  sensitive  to  the  degree  of  the  beam 
velocity  spread,  to  the  external  magnetic  field  profile,  and 
the  various  beam  and  circuit  parameters;  the  expression  for 
the  efficiency  was  derived  under  a  quasi-linear  approxima¬ 
tion.  Since  these  proposals  were  made  several  years  ago,  no 
experimental  results  have  yet  been  reported. 

Another  concept  for  increasing  the  bandwidth  of  the 
gyro-TWT  is  to  slow  down  the  wave  such  as  to  match  the  beam 
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axial  velocity  over  a  relatively  large  frequency  range .  In 
other'  words,  the  idea  suggested  here  is  to  find  a  structure 
capable  of  maintaining  the  grazing  condition  over  a  large 
enough  frequency  range.  Several  slow  wave  structures  have 
been  proposed,  Including  dielectric  loading  of  the  waveguide 
[24],  helix  loaded  waveguide  [25]-[26],  and  disc-loaded 
waveguide  £26]  —  [27]  .  No  experimental  results  on  these  ap¬ 
proaches  have  been  reported  either. 

Choe  and  Uhm,  [26],  investigated  the  disc  loaded  wave¬ 
guide  in  connection  with  the  gyrotron  application.  Using 
the  TEQ1  mode,  they  found  that  the  group  velocity  of  the 
wave  can  be  adjusted  by  varying  the  disc  inner  and  outer 
radii  a  and  b,  and  the  spacial  period  L  (see  figure  1.2  for 
notation) .  However,  the  TE01  mode  is  not  the  only  propagat¬ 
ing  mode  in  this  structure  and  hence,  mode  competition, 
[28],  is  a  possible  shortcoming  that  has  to  be  addressed.  On 
the  other  hand,  the  TEqj  mode  is  the  only  one  which  excites 
pure  TE  modes  at  the  discontinuities,  thus  simplifying  the 
analysis  of  the  fields  in  this  structure.  The  only  calcula¬ 
tion  given  in  [26]  on  the  gain  versus  frequency  showed 
finite  but  strongly  frequency  dependent  gain  over  a  larger 
than  usual  bandwidth  (13*).  It  is  our  understanding  that 
this  idea  is  not  pursued  at  NSWC  any  longer. 

In  the  present  effort,  we  also  analyse  the  periodical¬ 
ly  disc-loaded  cylindrical  waveguide  for  gyro-TWT  applica- 
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tion,  but  using  the  hybrid  mode.  The  principal  goal  of 

this  study  is  not  to  achieve  a  MM  wave  device,  but  to  achie¬ 
ve  a  device  for  S  to  K  band  range  applications  of  suffi¬ 
ciently  high  quality  that  it  can  compete  with  conventional 
TWT '  s  in  both  performance  and  cost.  The  choice  of  the 
hybrid  mode  shall  be  elaborated  further  later  on. 

Chapter  2  describes  the  general  approach  to  the  prob¬ 
lem.  It  describes  the  reasoning  leading  to  the  desired 
circuit  dispersion  characteristic.  In  chapter  3 ,  we  show 
that  the  cold  circuit  dispersion  relation  can  be  accurately 
described  by  a  polynomial  of  the  form 

ko2  =  a4k*  +  a2k^  +  ao 

in  the  first  Brillouin  zone.  We  make  use  of  the  polynomial 
approximation  for  the  cicuit  dispersion  relation  to  obtain 
the  gain,  small  signal  bandwidth  and  the  phase  characteris¬ 
tics  of  the  gyrotron  effect.  These  results  from  the  small 
signal  (linear)  analyses  were  obtained  in  the  two  cases  of  a 
spacial  growth  (complex  wavenumber,  real  frequency)  and  of  a 
temporal  instability  (complex  frequency,  real  wavenumber). 

In  chapter  4,  the  periodically  disc  loaded  cylindrical 
waveguide  is  analyzed.  The  dispersion  relation  obtained 
numerically  is  compared  to  the  experimental  results  whenever 
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possible.  The  results  are  then  fitted  with  the  polynomial 
approximation.  The  circuit  analysis  was  for  an  exp(j</>) 
azimuthal  variation  for  which  the  Floquet  modes  are  hybrid. 

Chapter  5  is  a  presentation  of  the  non-linear  analysis 
based  on  the  phase  averaging  technique  developed  by  Kuo  and 
Cheo,  £29]— [30] ,  and  adapted  to  the  present  structure. 
Saturated  gain,  bandwidth,  output  power  and  efficiency  are 
presented . 

The  last  chapter  compares  the  results  obtained  to  those 
published  for  different  configurations,  identifies  and 
discusses  some  pertinent  issues  for  further  work:. 

Two  appendices  are  included  for  convenience  of  the 
readers : 

*  Appendix  A  details  the  linear  dispersion  rela¬ 
tion  derivation  for  TE  modes  [31]  and  for  TM 
modes ; 

*  Appendix  B  is  the  derivation  of  the  non  linear 
int eg ro -differential  equation . 

A  list  of  references  is  provided  at  the  end. 
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CHAPTER  2 
GENERAL  APPROACH 

Most  of  the  works  on  gyro-TWT's  done  in  the  past  used 
structures  of  the  smooth  wall  type,  i.e.,  parallel  plate 
waveguides,  rectangular  or  cylindrical  waveguides.  All  these 
structures  have  in  common  the  fact  that  they  support  super- 
luminous  waves.  The  dispersion  relation  for  these  structures 
is  the  typical  hyperbola 

a >  2 

(-)  -  (k2  +  kc2)  =  0  (2.1) 

c 

where  kc  is  the  cutoff  wave  number  determined  only  by  the 
waveguide  geometry,  and  is  therefore  independant  of  both  the 
frequency  a>  and  the  propagation  constant  k. 

A  beam  propagating  in  a  helical  trajectory  along  the 
waveguide  will  interact  with  this  superluminous  wave.  The 
electrons  will  experience  forces  that  vary  as 

j(wt  -  kz) 

F  =  F0e 

If  the  beam  is  periodic,  and  its  velocity  in  the 
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transverse  plane  varies  as 

vt  =  vQcoskbz 

then  the  rate  of  change  of  the  energy  of  the  electrons  will 
be 


dE  j(o>t  -  kz) 

”  *  voFocoskbze 

dt 


=  "  v0F0 


j(«t  -  kz  -  kbz)  j(o>t  -  kz  +  kbz) 


+  e 


The  energy  exchange  shows  two  components,  but  only  the 
one  traveling  slower  than  the  wave  is  of  interest.  Indeed, 
when  the  slower  wave  is  traveling  at  about  the  same  speed  as 
the  electrons  in  the  longitudinal  direction,  cumulative 
bunching  and  energy  extraction  can  occur.  This  is  described 
by 


»t  -  (k  +  kb)z  =  (a>  -  kvz  -  n)t 

where  v2  *  z/t  is  the  electron  velocity  in  the  z-direct ion , 
and  kb  =  Q/vz  describes  the  periodicity  of  the  beam. 

Therefore,  when  the  Doppler  shifted  wave  frequency,  u>- 
kv^,  is  equal  to  the  cyclotron  frequency  0,  the  slower  wave 
of  energy  exchange  is  in  synchronism  with  the  beam:  this  is 
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known  as  the  grazing  condition.  This  description  does  not 
explain  the  need  of  actually  having  a  slightly  higher 
operating  frequency  than  the  Doppler  shifted  cyclotron 
frequency,  a  concept  known  in  the  community  as  A (d  >  0. 
However,  the  results  of  computing  the  hot  tube  dispersion 
relation  self-consistent ly  shall  indeed  show  that  the  modes 
with  gain  must  have  also  higher  operating  frequency, 
fulfilling  Aco  >  0. 

The  problem  is  that  this  synchronism  cannot  be  main¬ 
tained  over  a  large  frequency  band:  the  rapid  increase  of 
the  group  velocity  as  a  function  of  frequency  in  a  smooth 
waveguide  circuit  allows  only  a  very  narrow  band  interac¬ 
tion.  Therefore,  for  the  interaction  to  be  maintained  over  a 
larger  band,  the  group  velocity  should  not  vary  much  over 
all  the  required  band.  In  fact,  the  group  velocity  of  the 
wave  should  be  close  to  the  beam  velocity  in  the  longitudi¬ 
nal  direction  in  order  to  maintain  large  band  interaction. 
This  can  be  seen  through  the  following: 

the  beam  mode  is  defined  by 

a>  =*  kv2  +  U  (2.2) 

and  the  goup  velocity  of  the  wave  by  dro/dk  which  must 
be  equal  to  vz  to  maintain  synchronism.  That  is,  the 
best  possible  case  is  to  have  the  beam  z-velocity  and 


the  wave  group  velocity  equal  over  as  large  a  frequency 
band  as  possible. 

How  can  this  be  achieved? 

In  a  smooth  wall  waveguide,  the  transverse  wavenumber 
kc  is  a  constant,  determined  only  by  the  geometry  of  the 
guide,  and  not  by  either  m  or  k.  As  such,  it  is  not  possible 
to  match  the  waveguide  dispersion  relation  ( eq  2.1)  to  that 
of  the  beam  ( eq  2.2).  However,  if  we  can  make  the  transverse 
wave  number  a  or  k  dependant,  then  we  give  ourselves  the 
possibility  of  having  more  control  over  the  variations  of 
the  group  velocity  of  the  wave. 

In  the  derivation  of  the  smooth  wall  waveguide  disper¬ 
sion  relation  ( eq .  2.1),  the  left  hand  side  represents  the 
Helmholtz  operator  v2  +  (co/c)2.  If  an  artificial  boundary 
can  be  created  such  that  kc  is  no  longer  a  constant  but  <u  or 
k  dependant,  the  Helmholtz  operator  representation  on  the 
LHS  of  (eq  2,1)  is  replaced  by  a  different  function  f(w,k); 
and  f(w,k)  =  0  is  the  circuit  dispersion  relation.  Moreover, 
if  this  relation  can  be  matched  to  the  beam  dispersion 
relation  ( eq  2.2)  over  a  wide  frequency  range,  the  operating 
bandwidth  should  then  increase  accordingly.  We  need  there¬ 
fore  to  find  a  slow  wave  structure  so  that  the  group  velo¬ 
city  of  the  wave  matches  the  beam  velocity  down  the  wavegui¬ 
de.  A  circuit  that  can  achieve  the  two  stated  purposes 
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(  i .  e .  ,  a  transverse  wave  number  a>  or  k  dependant,  and  a 
group  velocity  about  constant  over  a  large  operating  band) 
may  be  found  in  the  form  of  a  periodic  structure. 

Indeed,  the  dispersion  relation  of  a  periodic  structure 
shows  (figures  2.1  and  2.2)  the  existence  of  a  region  in 
the  passband  where  w  versus  k  is  nearly  a  straight  line  (we 
are  interested  only  in  the  lowest  pass  band  and  in  the  first 
Brillouin  zone  in  this  region).  Choosing  the  beam  parame¬ 
ters  accordingly  (that  is,  in  such  a  way  as  to  make  the  two 
slopes  almost  equal  over  the  range)  should  solve  our  band¬ 
width  problem:  gain  can  be  achieved  over  a  larger  frequency 
band . 


Before  proceeding  further,  two  important  points  need  to 
be  made  : 

1-  In  a  periodic  structure,  the  rf  field,  as  per 
Floquet  theorem,  is  given  in  general  by 


00 

E  =  I  Ep ( r , 0 ) exp ( - j  kpz )  (2.3) 

P=  —  °° 

where 

kp  =  k  +  ( 2pn /L ) 

L  is  the  spacial  period 


Representing  the  Helmholtz  operator  on  the  LHS  of  (2.1)  by  a 
simple  algebraic  function  f(ai,k)  can  be  applied  to  the 
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fundamental  Floquet  mode  (the  term  p  =  0)  only.  It  is  shown 
ir>  chapter  4  that,  by  using  the  exp(j0)  mode  and  properly 
choosing  the  dimensions  of  the  structure,  it  is  possible  to 
have  the  beam  interact  strongly  only  with  the  p  =  0  term, 
and  minimally  with  the  other  terms  which  are  evanescent, 
decreasing  exponent ial ly  toward  the  axis  from  r  =  a.  Further 
more,  these  other  modes,  propagating  in  both  directions,  all 
have  much  slower  phase  velocity,  and  hence  do  not  interact 
effectively  with  the  beam. 

2 .  In  addition  to  TE  waves,  the  periodic  discon¬ 
tinuities  shall  also  introduce  TM  terms  which  should  in 
principle  be  taken  into  account.  That  is,  in  Ep  ( r ,  0 )  of 
equation  (2.3),  there  is  a  TM  part  as  well  as  a  TE  part.  For 
the  p  =  0  term,  the  combined  TE  and  TM  parts  form  the  modal 
field  which  interact  with  the  beam.  Thus,  contained  in  the 
linear  gyrotron  dispersion  relation,  should  also  be  the 
contributions  by  the  TM  part.  However,  the  TM  contribution 
is  negligible  and  becomes  zero  when  the  grazing  condition  is 
met*  . 


* 


S.P.Kuo,  private  communication. 
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CHAPTER  3 

SMALL  SIGNAL  ANALYSIS 


3.1  The  polynomial  approximation: 

In  the  previous  chapters  we  explained  how  it  is  pos¬ 
sible  to  overcome  the  limitation  on  the  bandwidth  of  a 
gyro-TWT  by  using  a  periodic  structure  instead  of  a  smooth 
walled  structure.  The  propagating  mode  in  this  circuit  will 
be  the  fundamental  TE^  mode  since  we  are  interested  in  the 
C/S  frequency  band.  However,  the  periodically  disc-loaded 
cylindrical  waveguide  which  we  will  investigate  for  gyro-TWT 
application  will  excite  all  of  the  Floquet  terms  of  the  TEln 
and  TMln  modes  for  an  exp(j0)  azimuthal  variation. 

In  such  a  circuit,  the  higher  Floquet  modes  (that  is 
for  p  ?  0)  do  not  interact  with  the  beam,  when  the  latter 
is  confined  near  the  axis  and  the  higher  modes  are  evanes¬ 
cent,  and  when  the  choice  of  the  waveguide  dimensions  are 
properly  made.  Representation  of  the  Helmholtz  operator 

at  2  to  2 

?2  +  (-)  by  (-)  -  (k2  +  kt2) 
c  c 

is  legitimate  when  only  one  single  mode  is  considered.  Such 
a  representation  can  also  be  legitimate  in  this  case,  except 
that  the  k^2  is  to  be  replaced  by  a  function  of  k  given  by 
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the  circuit  dispersion  relation.  Over  the  range  (o>,k)  of 
interest,  i.e.  the  first  Brillouin  zone,  this  can  be  well 
approximated  by  a  polynomial.  This  polynomial  should  not 
only  fit  the  actual  dispersion  curve  of  that  particuliar 
circuit,  but  should  also  satisfy  two  basic  conditions: 

1.  It  must  have  the  symetry  of  the  system.  A  periodi¬ 
cally  loaded  cylindrical  waveguide  show  an  even  symetry; 
the  polynomial  must  therefore  be  an  even  polynomial 

2.  At  k  =  tt/L,  the  slope  da>/dk  must  be  zero. 

Thus , 

W  2 

(-)  -  (a4k4  +  a2k2  +  aQ) 

c 

should  be  a  reasonable  form  for  the  polynomial* 

The  above  discussion  led  us  to  a  dispersion  relation 
for  the  periodically  loaded  cold  tube  which  is  simple  enough 
to  be  used  in  the  analysis  of  the  hot  tube,  and  yet  theore¬ 
tically  sound  to  be  reasonably  accurate. 

The  dispersion  relation  for  the  cold  tube  is  thus  given 
by 

«  2 

ko2  =  (“)  =  a4k4  +  a2k2  +  ao  (3.1) 

c 

The  three  constants  leading  to  the  approximate  circuit 
dispersion  relationship  can  be  found  in  any  appropriate  way. 
Generally,  two  conditions  should  be  met: 
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1.  aQ  =  cuc2/c2,  where  u>  c  is  the  lower  cut-off 
frequency,  obtained  either  analytically  or  experimentally. 

2.  do>/dk  =  0  at  k  =  n/L. 

The  third  relation  to  determine  all  three  coefficients 
can  be  set  by  some  curve  fitting  to  either  theoretical 
results  by  solving  the  boundary  conditions,  or  to  an  ex¬ 
perimental  curve.  If  neither  is  available,  any  reasonable 
assumption  may  be  tried,  e.g.,  d2w/dk2  =  0  at  k  =  7t/2L. 

The  polynomial  approximation  used  in  this  study  is 
shown  in  figure  3.4,  along  with  the  dispersion  curve  of  a 
smooth  walled  waveguide  of  the  same  outer  radius  b. 

However,  in  general,  the  polynomial  approximation  is 
found  from  the  numerical  results  of  the  complete  dispersion 
relation  described  in  chapter  4.  Experimental  plots  for 
different  values  of  a,  b,  and  L  are  presented  at  the  end  of 
chapter  4,  figures  4.1  and  4.2. 

3.2  Linear  analysis  of  the  hot  tube: 

In  the  litterature,  the  arrangement  of  the  gyro-TWT 
consists  of  an  annular  electron  beam  propagating  inside  the 
waveguide.  This  arrangement  is  appropriate  when  the  higher 
order  modes  are  the  ones  for  which  amplification  is  sought 
(i.e.,  the  interaction  of  interest  is  with  the  higher 
harmonics  of  the  electron  cyclotron  frequency) .  Indeed,  the 
beam  must  be  positioned  in  such  a  way  as  to  interact  with 
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the  rf  field  in  the  region  where  the  field  is  maximum,  a 
very  delicate  arrangement  in  some  cases:  figure  3.1  shows 
the  distribution  of  the  field  lines  for  a  waveguide  in  the 
whispering  gallery  mode.  However,  in  the  case  at  hand,  the 
fundamental  TE1;1  mode,  the  field  is  maximum  at  the  center  of 
the  guide  (figure  3.2);  hence,  the  beam  should  be  positioned 
at  the  center.  In  fact,  the  positioning  of  the  beam  is  not 
as  critical  as  when  considering  higher  order  interactions: 
appendices  A  and  B  show  that  the  energy  density  profile  in 
the  cross  section  of  the  guide,  the  geometrical  factor 
J ^ 2 ( 1 . 84RQ/RW )  changes  by  16%  only  when  RQ  is  changed  from  0 
to  30%  Rw .  Therefore,  an  electron  beam  located  near  the 
center  of  the  waveguide  with  beamlets  of  radius  RL  (Larmor 
radius)  centered  at  RQ  (the  guiding  center  radius)  equal  to 
RL  should  give  results  not  much  different  from  those  whose 
guiding  center  is  the  center  of  the  waveguide  (actually, 
Rq/rw  on  the  or^er  of  15%).  Thus,  the  arrangement  of  the 
gyro-TWT  in  the  fundamental  TE^  mode  is  depicted  in  figure 
3.3. 


The  electrons,  guided  by  a  uniform  magnetic  field  BQzQt 
move  along  helical  trajectories.  In  our  model,  we  assume 
that  the  beam  is  sufficiently  tenuous  that  its  space  charge 
electric  field  can  be  neglected  and  the  spacial  structure  of 
the  waveguide  mode  is  unaffected  by  the  presence  of  the 
beam,  the  beam  interacts  with  the  single  TE^  waveguide; 
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mode  (appendix  A  shows  that  the  additional  terms  due  to 
interaction  with  the  other  TE  and  with  the  TM  modes  can  be 
disregarded) ,  and  the  cold  circuit  dispersion  relation  used 
in  the  following  analysis  is  the  polynomial  approximation 
derived  previously  (eq.  3.1). 


For  the  disc-loaded  circuit,  the  TE  part  of  the  p  =  0 
term  has  the  same  field  distribution  as  that  of  the  TEj  ^ 
mode  of  a  smooth  walled  waveguide:  the  terms  on  the  RHS  can 
therefore  be  kept  as  such. 

The  small  signal  (linear)  dispersion  relation  for  the 
hot  tube,  with  the  above  assumptions  (see  also  the  discus¬ 
sion  in  appendix  A)  is  then: 
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1 

-j12(x)j1'(y)J2’(y)] 

2 

x  =  ktR0 

Y  =  ktRL 

1 

K11  =  <1 - 2>Jl2(ktRw) 

1 . 842 

k^  =  |^(w/c)2  -  k2j1//2 

i 

i 

At  this  point,  one  more  question  needs  to  be  answered. 
The  quantity  Rw  that  appears  in  equation  3.2  enters  into 
the  derivation  from  the  averaging  of  the  beam  current  over 
the  cross  section  of  the  smooth  waveguide.  It  is  the  boun¬ 
dary  where  or  1 (k^r)  is  zero.  Therefore,  the  periodi¬ 
cally  disc  loaded  cylindrical  waveguide  may  be  represented 
by  a  smooth  wall  cylindrical  waveguide  with  an  artificial 
conducting  boundary  at  R^  However,  this  boundary  is  not 
stationary,  but  function  of  frequency,  with  k^f^  =  1.84. 

Using  this  value  of  on  the  RHS  should  be  a  good  approxi¬ 
mation.  Since  the  per  unit  length  gain  of  the  tube  depends 
weakly  (1/3  power)  on  the  multiplying  factor  on  the  RHS  of 
equation  3.2,  the  error  introduced  should  therefore  be  very 
minor . 

The  first  term  on  the  right  hand  side  of  equation  3.2 
gives  the  gain,  and  the  second  term,  known  as  the  Weibel 
term,  is  a  damping  term.  Both  mechanisms  are  present  simul- 


taneously  and  compete  rather  than  reinforce  each  other.  The 
first  term  (the  one  multiplied  by  H^)  drives  the  instabi¬ 
lity  whereas  the  second  (the  one  multiplied  by  Qn)  has  a 
stabilizing  effect.  (Note:  when  =  0,  we  get  a  cubic 
equation  in  with  3  stable  roots).  The  Q  term  leads  to  a 
threshhold  beam  energy  beyond  which  the  CMI  can  be  excited: 

0t(°)  >  (critical ) 
i.e. ,  rt(0)  >  rt < critical ) 

We  see  that,  in  order  to  have  gain,  the  RHS  (which 
represents  the  coupling)  must  be  negative.  Complex  roots, 
whenever  they  occur,  give  the  gain  of  the  device.  Equation 
3.2  may  be  solved  for  either  w  or  k.  In  the  first  case,  k 
is  a  real  parameter  and  we  solve  for  complex  w,  w  =  wr  -  jei^ 
;  whereas  in  the  latter  case,  the  solution  is  a  complex 
value  for  k  =  kr  +  jk^,  a>  being  the  real  parameter.  Both 
approaches  were  used.  Graphs  1  through  8  are  plots  of  the 
results  of  the  linear  analysis.  It  needs  to  be  pointed  out 
that  the  polynomial  representation  of  the  cold  circuit 
dispersion  is  a  mathematical  approximation  in  the  first 
Brillouin  zone  of  the  true  dispersion  which  gives  all  bran¬ 
ches  in  the  entire  kQ  -  k  plane.  Thus,  it  does  not  contain 
all  the  physics  of  the  system.  Therefore,  in  solving  equa¬ 
tion  3.2  numerically,  if  any  solution  found  lies  near  the 
cold  circuit  characteristic  in  the  first  Brillouin  zone,  it 


it  can  then  be  regarded  is  legitimate.  Any  solution  found 
away  from  this  neighborhood  must  then  be  discarded. 

3«3  Discussion  of  the  Numerical  Results: 

In  the  following,  the  polynomial  approximation  used  for 
the  dispersion  relation  of  the  cold  tube  represents  the 
experimental  model  #8  os  described  in  chapter  4.  The  reason 
for  using  this  model  is  twofold: 

1.  The  dimensions  a,  b,  and  L  are  such  that  the  slope 
of  the  dispersion  curve  (which  is  equivalent  to  vz  = 
<3zc)  over  most  of  the  passband  is  small  enough  to  have 
positive  gain;  in  other  words,  the  energy  of  the  beam 
(i.e.,  r)  is  above  threshold. 

2.  With  reference  to  figures  4.1  and  4.2,  one  can 
readily  understand  the  compromise  one  has  to  make 
between  gain  and  bandwidth.  The  ultimate  bandwidth  is 
limited  by  the  difference  between  the  lower  (at  k  =  0), 
and  the  upper  (at  k  =  k/L)  cut  off  frequencies.  The 
larger  is  this  difference,  the  larger  the  bandwidth 
will  be.  But  this  also  leads  to  a  larger  group 
velocity:  less  time  is  available  for  the  wave  to  grow 
for  a  fixed  tube  length.  Further  more,  the  larger  the 
beam  parallel  velocity  leads  to  less  perpendicular  beam 
energy  available  for  the  interaction.  In  the 
computations  done  here,  a  fix  beam  voltage  of  60kV  is 
used.  This  implies  that  the  initial  r  = 


1  +  60/511 
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1.1174  and  (3Q  =  0.446.  The  most  recent  experimental 

results,  models  #6  through  #8,  showed  that  a  very  large 
range  of  dispersion  can  be  obtained,  from  pz  ~  0.08  for 
model  #6  (a  «  1.016cm,  b  =  2.93cm,  L  =  2.7cm)  at  the 
most  dispersive  end,  to  |3Z  «  1.0  for  the  earlier  models 
(a  =  2.22cm,  b  =  2.98cm,  L  -  3.3cm,  for  model  #2), 

almost  the  same  as  the  smooth  wave  guide.  A  desirable 
division  of  perpendicular  and  parallel  velocities  is 
about  vt/c  =  0.3  and  v2/c  =  0.33.  Model  #8  has  vg/c  - 
0.32. 

In  this  analysis,  we  therefore  used  the  values  of 
experimental  model  #8: 

L  =  0.027m 

b  =  0.0298m 

a  =  0.5b  =  0.01486m 

With  reference  to  figure  3.4,  note  how  remarkable  is 
the  agreement  between  the  experimental  data  and  the 
polynomial  approximation  used.  This  proves  our  point:  the 
dispersion  relation  of  this  structure  represented  by  an 
infinite  by  infinite  determinant  can  be  accurately  described 
by  equation  3.1. 

The  coeficients  in  equation  3.1  are: 
a4  =  -2.74X10"5 
a2  =  0.7413 

a0  = 


8046.1 
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The  initial  parameters  for  the  beam  were  the  following 
VB  =  60kV 
IB  =  5A 
r  =  1.1174 

Ro  =  RL 
Bq  -=  1500G 

3,3-3  numerical  Results: 

Several  beam  parameters  were  used  in  order  to  compare 
the  performances  of  their  respective  tubes  (  Fig  .  3 . 5  )  .  Two 
beams  are  then  used  for  further  linear  and  non  linear 
studies:  beam  B1  which  meets  the  grazing  condition  over  the 
largest  range  of  frequencies,  and  beam  B3  which  achieves  the 
flatest  linear  gain  over  the  passband ( Fig . 3 . 6 ) . 

Figure  3.5  shows  the  cold  circuit  dispersion  curve  (WG) 
along  with  three  straight  lines  that  describe  the  dispersion 
curves  of  the  three  beams  used  (B1-B3).  These  three  beams 
differ  by  the  value  of  the  electron  cyclotron  frequency  Q, 
and/or  their  slope  AkQ/Ak;  the  beam  voltage  and  the  beam 
current  are  kept  the  same  for  all  three  beams,  i.e.,  Vb  = 
60kV  and  Ib  =  5A. 

Figure  3.6  is  a  plot  of  the  gain  versus  frequency  for 
all  three  interactions  WG~B^  using  the  complex  a> ,  real  k 
computer  program.  The  WG-B1  interaction  shows  the  largest 
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small  signal  gain  of  all  three  interactions. 

We  are  not  however  interested  in  getting  the  largest 
gain  because  it  may  give  rise  to  some  absolute  instability. 
What  is  more  desirable  is  a  reasonable  high  and  uniform  gain 
over  a  large  frequency  range,  that  is,  we  would  like  to 
obtain  a  gain  as  flat  as  possible  over  as  large  a  frequency 
range  as  possible.  The  reason  for  this  is  to  have  the  tube 
output  the  maximum  (or  close  to  the  maximum)  power  over  that 
range  of  frequencies,  hence,  to  have  a  tube  working  at 
maximum  efficiency  over  a  relatively  large  bandwidth.  The 
plots  in  figure  3.6  point  obviously  to  B3  to  achieve  the 
best  performance. 

As  mentioned  earlier,  equation  3.2  may  be  solved  by 
either 

1.  Let  k  be  a  real  parameter,  then  looking  for 
roots  of  a) ;  or 

2.  By  letting  a>  be  a  real  parameter  while  looking 
for  roots  of  k. 

Whenever  complex  roots  occur,  gain  is  indicated.  In  the 
first  approach,  an  equivalent  per  unit  length  gain  (k^)  or 
phase  shift  (kr)  may  be  obtained  by  letting  k^  =  «jj/vg  and 
kr  =  o»r/Vg.  The  equation  is  4th  degree  in  w.  It  is  found 
that  over  a  large  range  of  k,  two  complex  roots  can  be 
obtained.  Two  other  roots  are  real,  and  hence  are  stable 
modes.  The  two  complex  roots  are  conjugate  to  each  other. 
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modes.  The  two  complex  roots  are  conjugate  to  each  other. 
Figures  3 . 7 ( f or  the  WG-B1  interaction)  and  3. 12 (for  the  WG- 
B3  inteaction)  show  the  kr  solutions  versus  k  as  compared 
with  cold  circuit  characteristics  and  the  beam  characteris¬ 
tic.  It  is  of  interest  to  note  that: 

1.  The  kr  of  the  complex  mode  lies  above  the  beam 
characteristic,  and  is  almost  parallel  to  it.  This  shows 
that  Aw  >  0 ,  and  the  group  velocity  is  nearly  constant  and 
equal  to  vz  of  the  beam. 

2.  The  two  stable  modes  are  below  the  beam  charac¬ 
teristic  ( Aw  <  0).  The  one  immediately  below  is  a  modified 
waveguide  mode  with  a  smaller  phase  velocity.  The  other 
stable  mode  is  a  backward  wave;  it  is  however  too  far  from 
the  first  Brillouin  zone  to  be  taken  seriously. 

3.  Figure  3.12  shows  that  for  low  frequencies,  no 
instability  is  excited:  the  solution  of  equation  3.2  shows 
only  stable  roots.  Then,  at  some  frequency,  two  of  the  real 
roots  merge  into  two  complex  conjugates  whose  real  part 
presents  the  characteristic  feature  of  the  instability  been 
excited,  i.e..  Aw  >  0. 

Figures  3.8  and  3.9  (for  the  WG-B1  interaction)  and 
figures  3.13  and  3.14  (for  the  WG-B3  interaction)  show  the 
gain  (in  dB/m)  versus  frequency  obtained  by  the  two  approa¬ 
ches.  Figure  3.10  and  3.15  compares  the  two  sets  of  results: 
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the  agreement  is  remarkable,  especially  for  the  WG-B3  inter¬ 
action. 

Figures  3.11  and  3.16  show  the  per  unit  length  phase 
delay  versus  frequency,  showing  the  linear  phase  relation  so 
very  important  for  amplifier  applications. 

To  summarize,  the  linear  analysis  of  the  gyro-TWT  shows 
that  reasonably  high  and  uniform  gain,  as  well  as  an  almost 
linear  phase  shift  can  be  achieved  with  this  design.  It  also 
predicts  a  20.4%  bandwidth. 


Positioning  of  beam  modes  relative  to  the  waveguide  mode 
(WG)  represented  by  the  polynomial  approximation 


4. 
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5. 


5.5  (GHz) 


Fig.  3.10 

Figures  3.8  and  3.9  superimposed  (WG-B1  interaction) 


WG-B1  interaction 
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CHAPTER  4 

THE  PERIODICALLY  DISC  LOADED 
CYLINDRICAL  WAVEGUIDE 


4.1  Derivation  of  the  dispersion  relation: 

The  disc  loaded  waveguide  circuit  has  been  analyzed  in 
the  past  in  connection  with  studies  of  the  conventional  TWT 
[24];  only  azimuthally  symmetric  TM  modes  were  of  interest. 
Choe  and  Uhm  [23]  analyzed  the  axially  symetric  TE  modes  of 
the  structure.  These  are  the  only  cases  where  pure  TE  or 
pure  TM  modes  are  excited.  In  this  section,  we  solve  the 
boundary  value  problem  of  the  cold  disc  loaded  circuit,  thus 
leading  to  the  circuit  dispersion  relation. 

The  discs  in  our  structure  excite  all  TEln  and  all 
™ln  modes •  They  are  however  evanescent,  except  for  the 
dominant  TE^^  mode.  Hence,  for  a  given  choice  of  beam  para¬ 
meters,  only  the  dominant  mode  will  couple  to  the  beam, 
thus  amplified. 

The  expression  for  the  fields  will  therefore  consist  of 
all  the  Floquet  terms  of  both  TE  and  TM  modes.  With  respect 
to  figure  1.2,  and  using  the  exp(ja>t)  convention,  the  fields 
for  an  exp(j</>)  azimuthal  variation  can  be  expressed  as: 


For  region  I  (r  <  a) 
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For  region  II  (a  <  r  b) 
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Where  Xp  and  r?m  become  imaginary,  the  Bessel  functions 
are  replaced  by  the  modified  functions:  Jn(x)  ->  In(x), 

Yn(x)  ->  Kn(x) . 

We  define 

zE<*»Br>  =  Jl^mr> - T - Yl^mr) 

Y1  ^mb) 
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(4.3) 
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and  Ap,  Cp '  are  column  vectors  representing  the  amplitudes 
of  the  fields. 

For  non  trivial  solution  (Ap,  Cp r ) ,  the  determinant  of 
the  matrix  must  be  zero:  this  is  then  the  dispersion  rela¬ 
tion  of  the  cold  circuit  for  the  entire  Brillouin  diagram. 

We  would  like  to  emphasize  here  the  fact  that  the 
interaction  between  the  beam  and  the  hybrid  wave  is 
equivalent  to  the  interaction  between  the  beam  and  a  TE^ 
smooth  walled  waveguide.  Indeed,  equation  4.3  shows  that  a 
beam  located  near  the  axis  will  interact  mostly  with  the  p  = 
0  term  of  the  wave  fields,  the  other  terms  having  but  a 
small  effect  on  the  fields  near  the  axis.  This  is  so  because 
the  Bessel  functions  will  have  real  argument  only  for  the  p 
=  0  term,  and  imaginary  argument  for  all  the  other  terms; 
meaning  that  the  field  distribution  near  the  axis  will  look 
very  much  like  the  one  for  the  TE^  mode  since  the  modified 
Bessel  function  goes  to  0  towards  the  axis.  Another  reason 
is  that  the  phase  velocity  of  the  p  ^  0  waves  is  much  less 
than  the  beam  velocity,  hence,  even  if  some  interaction 
takes  place,  it  cannot  sustain  an  instability.  The  hybrid 

field  distribution  at  the  beam  location  is  therefore 
about  the  same  as  that  of  a  TE^  smooth  walled  waveguide 
mode,  a  fact  that  allows  us  to  use  the  RHS  of  the  linear 
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dispersion  relation  (eq.  3.2)  as  such,  adjusted  for  this 
case  (i.e.,  substituting  k^2  by  kD2  -  k2). 

Of  course,  this  is  not  the  only  possible  approach. 
Other  approaches  are  possible,  and  indeed  are  used.  They 
however  have  the  disadvantage  that,  although  they  give 
correct  results,  the  physical  picture  is  not  self-evident, 
imbedded  too  deeply  in  the  mathematics. 

4.2  Discussion: 

The  dispersion  relation  in  this  case  is  very  involved: 
it  is  the  solution  of  a  double  infinite  determinant  erual  to 
zero,  which  is  of  course  impossible  to  compute.  Truncating 
this  determinant  to  some  manageable  size  will  undoubtedly 
introduce  some  error  which  has  to  be  minimizeid.  The  jus¬ 
tification,  if  at  all,  for  truncation  is  that  the  modes 
above  a  certain  level  are  of  minimal  influence  to  the  resul¬ 
ts.  However,  .by  doing  so,  it  actually  means  that  the  modes 
discarded  are  of  no  importance,  while  those  kept  are  of 
equal  importance.  Unfortunately,  neither  of  this  is  true. 
Not  all  lower  modes  are  of  equal  importance  to  the  dominant 
first  TE  mode.  Meanwhile,  the  sharp  edges  at  the  discs  are 
expected  to  cause  slow  convergence  at  the  higher  modes.  On 
the  other  hand,  very  often,  using  only  one  mode  is  suffi¬ 
cient  to  give  reasonable  results,  such  as  in  magnetron 
cavity  analysis.  Which  means  that  the  determinant  should  not 
be  truncated  by  just  eliminating  the  higher  order  term. 
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This  and  other  similar  problems  of  truncation  have  faced  the 
electromagnetic  community  for  decades.  Numerous  methods  for 
dealing  with  this  problem  have  been  devised,  none  of  them 
however  giving  satisfactory  results  for  all  cases.  In  this 
analysis,  very  often  an  error  of  10%  compared  to  experiment 
is  regarded  as  satisfactory.  Whereas  in  this  study,  whenever 
we  can  devise  an  approach,  the  error  is  of  the  o  der  of  1% . 
We  were  fortunate  however  that  a  parallel  effort  was  under¬ 
taken  in  which  the  dispersion  characteristic  of  disc-loaded 

* 

circuits  were  experimentally  measured  :  eight  models  were 
built  and  measurements  were  performed.  Figures  4.1  and  4.2 
are  the  plots  obtained. 

With  regard  to  the  numerical  approach  used,  since  we 
have  no  way  of  knowing  a  priori  which  terms  are  of  lesser 
influence,  a  trial  and  error  approach  was  used,  the  results 
compared  to  those  given  by  the  experiments.  As  we  mentioned 
earlier,  by  properly  choosing  the  value  of  the  axial  period 
L,  it  is  possible  to  have  the  beam  interact  strongly  with 
the  fundamental  (p  =  0)  Floquet  mode  only.  All  the  other 
forward  Floquet  terms  are  evanescent,  decreasing  exponen¬ 
tially  towards  the  axis  from  r  =  a.  The  choice  ^or  L  is 
therefore  dictated  by 


Experiments  were  performed  by  G.  Wurthmann,  J.  Heary 
and  C.  Bates  of  Electronic  Device  and  Technology  Laboratory 
(ETDL)  at  Fort  Monmouth.  The  results  were  communicated  to  us 
privately . 
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for  p  =  0 
for  p  /  0 


(4.5) 


where  kp  is  given  by  equation  4.2. 

For  a  fixed  range  of  frequencies,  the  value  for  L  can 
be  determined  from  equations  4.5  and  4.2.  For  the  frequency 
range  of  interest,  i.e.,  such  that  the  lower  cutoff  frequen¬ 
cy  3. GHz  <  fc  <  4 . GHz ,  it  is  found  that  L  <  .025  m. 

The  value  for  b  is  determined  by  the  lower  cut-off 
frequency  such  that  k^b  =  1.84,  kt  being  the  wave  number 
corresponding  to  a  3. GHz  smooth  wall  WG  cut-off  frequency. 

By  changing  the  size  of  the  discs,  the  slope  of  the 
dispersion  curve  can  be  adjusted,  as  can  be  seen  in  graphs 
9  and  10.  With  a  choice  of  kt  =  62.83,  the  value  for  b  will 
be  2 . 93cm . 


It  turns  out  that  the  value  for  L  is  also  of  paramount 
importance  to  the  numerical  analysis:  for  the  value  of  L 
such  that  equation  4.5  is  satisfied,  we  always  have  <  o. 
In  this  case,  whatever  we  do,  the  numerical  results  are 
erroneous:  we  get  a  complete  dispersion  curve  well  above 
where  it  should  be.  In  fact,  the  numerical  simulation  yields 
results  only  when  rj  L  is  real,  that  is,  when 


? 
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^l2  =  ko2  "  ('t/L)2  >  0 

If  we  choose  L  such  that  kQ  >  n/L  for  part  of  the  dispersion 
curve,  then  we  get  very  accurate  results  for  that  portion: 
indeed,  in  this  case,  the  piece  of  the  dispersion  curve 
above  kQ  =  nf  L  matches  almost  perfectly  the  experimental 
curve*.  For  the  piece  below  kQ  =  ;c/L,  the  results  are  howe¬ 
ver  far  off,  well  below  what  the  measurements  give  us.  So 
far,  we  do  not  have  any  satisfactory  explanation  for  this 
numerical  problem.  We  do  however  believe  that  we  are  not  far 
from  reaching  a  solution.  For  example,  Galerkin's  approach 
is  now  being  looked  into. 
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CHAPTER  5 


NON  LINEAR  ANALYSIS 


5 . 1  Background : 

The  linear  analysis  provides  us  with  a  wealth  of  infor¬ 
mation  with  regard  to  the  smail  signal  performance  of  the 
tube.  However,  according  to  the  linear  theory,  if  the  wave 
is  gaining  energy,  it  will  keep  on  gaining  energy  forever. 
Indeed,  there  is  no  provision  in  the  linear  theory  for  any 
kind  of  saturation  mechanism  that  the  wave  must  experience. 
In  reality,  some  mechanism  is  always  present  to  limit  the 
output.  To  assess  the  performance  of  a  power  amplifier,  a 
non  linear  theory  must  then  be  developed,  which  must  predict 
the  saturation  level,  as  well  as  all  the  other  important 
performance  results  such  as  the  saturation  length,  the  power 
output,  the  efficiency  of  the  tube. 

Two  different  but  simultaneously  present  (in  the  fun¬ 
damental  interaction)  mechanisms  are  responsible  for  the 
saturation  : 

-  Depletion  of  free  energy 

-  Loss  of  phase  synchronism  (phase  trapping) 

Depletion  dominates  when  the  beam  energy  is  slightly 

above  threshhold.  Saturation  occurs  as  soon  as  the  beam 
loses  a  small  amount  of  energy  and  the  system  becomes  line- 
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arly  stable.  This  effect  becomes  also  apparent  when  the 
electrons  have  lost  enough  energy  to  the  wave  so  that  is 
no  longer  high  enough  to  cause  amplification.  Phase  trapping 
occurs  when  the  beam  energy  is  well  above  threshold:  an 
average  electron  loses  so  much  energy  that  its  relativistic 
cyclotron  frequency  no  longer  matches  the  wave's  to  favor 
unstable  interaction . 

Both  mechanisms  are  important  for  the  fundamental 
cyclotron  harmonic  interaction.  However,  the  dominant  mecha¬ 
nism  depends  on  the  initial  choice  of  the  beam  parameters: 
the  maximum  free  energy  per  particle  available  to  the  wave 
is 


E  =  [r(0)  -  r ( critical ) j mQc2 

Non-linear  behavior  of  the  gyrotron  mechanism  has 
mostly  been  done  by  particle  simulation  [37]— [38] ,  and  a 
combination  of  both  analytic  and  simulation  in  part  [34]. 
These  approaches  need  large  computers  and  are  hence  expen¬ 
sive  if  large  amounts  of  data  is  needed  for  trial  and  error 
design  work.  In  the  present  effort,  we  adapt  the  phase 
averaging  technique  by  Kuo  and  Cheo  [29],  and  by  Kuo  et  al  . 
[30]  to  the  present  structure.  In  [29]  and  [30]  a  single 
third  order  ordinary  differential-integral  equation  was 
derived  for  the  slow  varying  amplitude.  The  approach  is  less 
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accurate  than  the  particle  simulation  but  readily  solvable 
on  a  personal  computer.  Thus  it  is  most  convenient  for 
estimating  the  ultimate  performance  of  a  power  amplifier. 
Appendix  B  gives  the  essence  of  the  efforts  in  [29]  and 
[30].  The  basic  equation  is  given  by: 
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t 

Is  *  J  dt 1  a(t')  sin<A0(t  -  tf)> 
o 

Note  that  linearizing  equation  5.1,  we  find  that  in 
order  to  have  gain,  we  need  to  have 

(u)2  -  k2c2)(3t2H11  (x,  y)  - 

(a)  -  kvz  -  n)(w  -  kvz)Q11(x,y)  >  0 

i.e.,  the  linear  condition. 


The  output  power  of  the  tube  was  derived  through  the 
integration  of  the  Poynting  vector  over  the  cross-section  of 
the  tube. 
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Using  the  following  relationship 


J  rJm2(ktr)dr  = 
o 

R  ^ 

~  Jm-1  (  kt**w )  Jm+1  (kt^w  ^  ] 
and  substituting,  we  finally  get: 
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5.2  Numerical  Results 

For  a  homogeneous  cylindrical  waveguide  (refer  to 
equation  5.1), 
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This  term,  which  comes  from 

/  rJm2(ktr)dr 
o 


is  a  geometrical  factor  that  represents  the  way  the  power 
density  is  dis  tr  ibuted  over  the  cross* sec t ion  .  For  the 
geometry  considered,  this  term  should  not  be  very  much 
different  from  the  one  of  the  cylindrical  guide  (see  chapter 
3  for  a  discussion  about  1^) . 

The  above  integral -di f ferent ial  equation  was  then 
integrated  using  the  solution  k  and  w  of  the  linear  problem. 
Figures  5.1  (for  the  WG-B1  interaction)  and  5.6  (for  the  WG- 
B3  interaction)  show  the  wave  amplitude  EQ  versus  the  inter¬ 
action  length  for  several  frequencies*.  The  wave  amplitude 
at  the  input  of  the  tube  was  the  same  in  both  cases,  and 
chosen  very  small  so  that  the  results  of  the  linear  analysis 
apply  at  the  tube  entrance.  Several  conclusions  can  be 
inferred  from  these  two  figures: 

1.  In  the  WG-B 1  case,  the  growth  rate  of  the  wave  is 
approximately  linear,  except  near  the  saturation  point, 
and  the  gain  per  unit  length  in  the  pass-band  is  about 
the  same  as  the  small  signal  gain  derived  in  the  linear 
analysis,  which  proves  the  consistency  of  this  .appro¬ 
ach.  This  also  means  that  the  linear  approach  can  bo 
used  even  for  large  values  of  the  wave  intensity,  as 
long  as  we  are  far  enough  from  the  saturation  level 
the  tube. 
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real  k  results  were  used 


of 


58 


The  WG -B3  case  is  different  in  the  sense  that  the 
growth  rate,  although  linear,  is  different  from  the  one 
obtained  in  the  linear  analysis.  The  reason  is  that  the 
beam  and  the  wave  are  far  apart,  and  therefore  the 
nonlinear  relationship  is  different  from  the  one  deri¬ 
ved  in  [29]  &  [30].  However,  when  the  two  dispersion 
relations  get  closer  to  each  other,  i.e.,  near  the 
higher  cut  off  frequency,  both  analyses  yield  the  same 
results . 

2.  The  saturation  value  of  the?  electric  field  intensity 
(i.e.,  the  value  of  EQ  at  saturation)  decreases  as  the 
frequency  is  increased.  Also,  the  saturated  electric 
field  intensity  decreases  monotonical ly  as  the  frequen¬ 
cy  is  increased,  then  sharply  near  the  higher  cut  off 
frequency . 

3.  The  higher  the  frequency,  the  longer  it  takes  the 
wave  to  reach  saturation. 

4.  The  maximum  value  of  the  saturated  electric  field 
intensity  always  occurs  at  the  lowest  frequency  * . 


* 

This  relates  to  the  experimentally  observed 
the  largest  linear  gain  is  obtained  when  the 
cylotron  frequency  is  closest  to  the  lower  cut  off 
of  the  guide. 


fact  that 
elect  ron 
frequency 
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Figure  5.2  (WG-B1  case)  is  a  plot  of  the  maximum  output 
power  available,  i.e.,  the  saturation  power,  versus  frequen¬ 
cy.  This  power  increases  as  frequency  is  increased  and  shows 
a  sharp  drop  near  the  higher  cut  off  frequency. 

Figure  5.3  show  the  saturation  efficiency  versus 
frequency.  We  defined  this  efficiency  as  the  ratio  of  the 
wave  power  to  the  dc  power  of  the  beam  (the  Vg  [g  product). 
The  saturation  efficiency 


used  by  many  authors  is  meaningless  in  this  case  because  it 
is  greater  than  1.  This  is  due  to  the  averaging  process  used 
to  derive  the  integral-differential  equation.  Taking  more 
terms  into  cons idera t ion  when  averaging  *  the  initial  phases 
of  the  electrons  would  certainly  be  more  accurate,  at  the 
expense  however  of  having  to  solve  a  5th  order  differential- 
integral  equation  at  least* 

For  a  complete  final  picture  of  the  gyrotron  amplifier, 
several  choices  have  to  be  made:  indeed,  a  trade-off  between 
bandwidth,  output  power  and  tube  length  is  required. 

For  example  (see  figure  5.1),  a  20dB  tube  at  4.497GHz 
"center  frequency"  (defined  as  the  one  with  a  saturated 
output  electric  field)  will  be  37cm  long;  for  the  same  input 
electric  field,  a  wave  at  another  frequency  would  have 
either  passed  the  saturation  point  or  would  still  be  gaining 
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energy  from  the  beam.  In  this  example,  all  the  frequent;  ies 
below  4.497GHz  will  be  at  the  saturation  point  and  the 
higher  frequencies  are  close  to  it.  This  therefore  repre¬ 
sents  the  best  operating  conditions  for  the  tube  since  the 
power  output  will  be  maximum  (or  close  to  the  maximum),  and 
we  also  get  a  substantial  24.9%  bandwidth  (see  figures  5.4 
and  5.5)  .  The  bandwidth  can  be  defined  as  the  range  of 
frequencies  such  that: 

1.  the  output  electric  field  EQ  >  E0max/V2 

2.  the  output  power  PQ  >  PQmax/2. 

In  this  case,  they  do  not  represent  the  same  thing. 

Another  example  would  be  to  choose  the  "center  frequen¬ 
cy"  of  the  2  OdB  tube  to  be  5.259GHz.  But  then,  the  tube 
would  be  42cm  long  with  most  of  the  lower  frequencies  past 
the  saturation  point,  hence  a  lower  power  output  and  a 
smaller  bandwidth.  This  situation  arises  if  one  chooses 
maximum  saturation  power  frequency  as  the  "center  frequency" 


of  the  tube. 


p 

<kW> 
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Saturation  Power  vs  frequency 
WG-B1  interaction 


Saturation  Efficiency  vs  frequency 
WG-B1  interaction 


Output  Power  vs  frequency  for  the  37cm  long  tube 


1.6MV/m  at  4.283GHz 
1.6MV/m  at  4.497GHz 
J5MV/m  at  5.5  GHz 


BWS  24.9% 


Fig.  5.5 

The  20dB  tube  characteristics  (WG-B1  case) 
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CHAPTER  6 
CONCLUSION 

The  results  of  the  linear  and  non- linear  analyses  show 
that  a  periodically  disc-loaded  cylindrical  waveguide  design 
can  be  a  good  contender  for  gyrotron  amplifier  applications. 

A  37cm  long  tube  (20dB  saturated  gain)  can  deliver  a 
peak  output  power  of  about  48kW  at  5 . 3GHz  with  a  18%  ef¬ 
ficiency.  A  nearly  straight  line  phase  characteris tic ,  a 
bandwidth  of  about  25%  and  the  fact  that  the  structure  is 
non-dispersive  make  it  suitable  for  wide  band  power  amplifi¬ 
cation  in  communication  /radar. 

The  only  tubes  that  have  been  experimentally  tested  in 
-the  same  frequency  band,  are  two  Varian  tubes,  [18*]-[19]. 
They  were  both  of  a  smooth  wall  cylindrical  design;  however, 
the  magnetic  field  profile  in  the  second  tube  was  not  uni¬ 
form,  varying  from  1800G  at  the  input  of  the  tube  to  1900G 
at  the  output,  and  also  some  distributed  losses  were  in¬ 
cluded  in  this  second  tube.  The  first  tube  operated  at  a 
beam  voltage  =  60kV  and  a  beam  current  1^  =  5A,  whereas 
the  second  operated  at  Vfe  =  65kV  and  Ife  =  8A.  Both  tubes 
were  43cm  long,  for  an  18dB  saturated  gain  for  the  first 
tube  and  20dB  saturated  gain  for  the  second.  The  small 
signal  gain  was  approximately  the  same  in  both  cases,  24dB/m 
for  the  first  tube,  26dB/m  for  the  second.  The  peak  power 
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t: 


output  achieved  was  50kW  at  5GHz  (with  16.6%  efficiency)  for 
the  first  tube,  and  128kW  at  5 . 2GHz  (with  24%  efficiency  for 
the  second  tube;  the  bandwidths  were  respectively  6%  and 
8,8%.  Phase  linearity  for  the  second  tube  was  within  —  26°; 
it  was  not  reported  for  the  first  tube. 

In  all  the  work  performed  for  this  study,  the  beam  was 
assumed  to  be  monoenergetic,  the  space  charge  effect  and  the 
velocity  spread  were  both  neglected.  Chen  and  Chu,  [44], 
show  that  the  gain  is  insensitive  to  the  velocity  spread  in 
the  longitudinal  direction  for  monoenergetic  electrons,  for 
Ap2/pzo  =  0,  25%,  50%;  also,  according  to  their  study,  in 
regions  where  the  growth  rate  is  small,  a  spread  in  vz  can 
enhance  the  growth  rate'.  These  results  contradict  previous 
ones  (Caplan,  Lin,  Chu,  1982  as  reported  by  Feinstein  and 
Felch,  1987,  [43])  stating  that  the  output  efficiency  is 
degraded  by  a  factor  of  2  from  that  achievable  with  an  ideal 
beam  when  Apz/pz0  =  15%. 

In  that  same  paper,  Chen  and  Chu  confirm  the  existence 
of  a  space  charge  instability,  and  show  that  it  is  an  elec¬ 
tro-static  cyclotron  instability  that  exists  in  the  beam, 
independently  of  the  circuit.  It  can  lead  to  significant 
noise  enhancement  in  gyrotrons  and  gyrotrons  amplifiers.  The 
noise  level  was  measured  by  Ferguson  et  al .  to  be  in  44~52dB 
range  above  thermal,  comparable  to  the  small  signal  gain. 


i 
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These  problems  have  not  been  addressed  in  this  work, 
which  was  primarily  concerned  with  the  feasability  and 
characteristics  of  a  high  quality  gyro-TWT  using  a  slow  wave 
structure  as  its  circuit.  The  upcoming  experiments  at  Fort 
Monmouth  with  a  hot  tube  based  on  this  design  should  yield  a 
wealth  of  reliable  information. 


NOTE  ABOUT  THE  APPENDICES 


In  appendix  A,  the  derivation  of  the  linear  dispersion 
relation  is  presented  in  detail.  The  result  has  been 
published  first  in  [37]  and  subsequently  used  by  several 
investigators.  The  derivation  has  never,  however,  been 
published.  It  has  been  included  here  for  the  benefit  of  the 
newcomer  to  this  field  and  for  the  sake  of  completeness . 

The  detailed  derivation  of  the  non-linear  equation 
describing  the  evolution  of  the  radiation  fields  has  been 
included  in  appendix  B  for  similar  reasons.  It  was  first 
derived  by  S.P.  Kuo  and  B.R.  Cheo  [29]. 


APPENDIX  A 


DERIVATION. OP  THE  LINEAR  DISPERSION  RELATION 

A.l  General: 

We  have  seen  that  electrons  subjected  to  rapidly  fluc¬ 
tuating  electromagnetic  field  will  deviate  from  their  unper¬ 
turbed  orbits.  The  new  motion  can  be  resolved  into  a  quasi¬ 
static  term  plus  a  rapidly  oscillating  term.  It  is  the 
slowly  varying  part  of  the  motion  which  is  responsible  for 
the  instability  and  subsequent  nonlinear  behavior. 

The  arrangement  of  the  gyro-TWT  consists  of  an  annular 
electron  beam  propagating  inside  the  disc  loaded  waveguide, 
The  electrons,  guided  by  a  uniform  magnetic  field  B0zQf  move 
along  helical  trajectories.  In  this  model,  we  assume  that 
the  beam  is  sufficiently  tenuous  so  that  its  space  charge 
electric  field  can  be  neglected  and  the  spacial  structure  of 
the  waveguide  mode  is  unaffected  by  the  presence  of  the 
beam.  The  beam  interacts  with  a  single  TEmn  waveguide  mode, 
where  m  and  n  are  respectively  the  azimuthal  and  radial 
eigen  mode  numbers.  The  dynamics  of  the  electron  beam  is 
described  by  the  Vlasov  equation. 

In  this  appendix,  we  start  with  the  derivation  of  the 
linear  dispersion  relation  of  the  gyrotron  first  for  a  beam 
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interacting  with  a  TEmn  mode,  then  for  a  beam  interacting 
with  a  TMmn  mode,  combining  at  the  end  the  two  for  a  general 
linear  dispersion  relation. 

It  is  then  shown  that  the  term  due  to  the  interaction 
beam-TMmn  is  much  smaller  than  the  one  due  to  the  beam-TEmn 
interaction  and  hence,  can  be  discarded  without  losing  much 
information  regarding  the  behavior  of  the  hot  tube.  Siraiia- 
rily,  the  terms  due  to  the  non-fundamental  interaction  beam- 
TEmn  have  a  much  smaller  contribution  than  the  one  due  to 
the  dominant  mode. 

A. 2  Derivation  of  the  linear  dispersion  relation: 

Let's  start  with  Vlasov  equation 

[  it  +  ?i-'  -  e(E  +  u  K  "  0 

where  f  is  the  electron  distribution  function,  and  let's 
define  in  the  linear  approximation  f  =  fQ  +  fj,  where  fQ  is 
the  initial  and  fj  the  perturbed  distribution  functions 
respectively,  with  <<  f  .  Then 

[  It  +  rm  ,V  '  e(V  *  Bozo)-vp  ]  fl  =  e<Ei  +  v  x  Bl> -vpfo(A,1) 

where  we  have  neglected  on  the  right  hand  side  the  second 
order  small  term  [e(Ej  +  v  x  B^)  Ei  an<^  Bi  represent 
respectively  the  wave  electric  field  and  magnetic  displace¬ 


ment  . 
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To  solve  Vlasov  equation,  we  use  the  Lagrangian  coor¬ 
dinate  system,  and  define  an  unperturbed  trajectory 


dr  p 

cat  rm0 


dp  p 

—  =  -  e  v  x  zQ  =  -  e - x  BQ  zQ 

dt  rmQ 


(A. 2) 


n  p  x  zQ 


where  0  =  eB0/rm0  is  the  relativistic  electron  cyclotron 

frequency. 

Then,  along  the  trajectory: 

d  p  a  dr  dp  d 

--  + - .v  -  e  ( v  x  Bq  zq)  =  --  +  v  +  — .  v  =  -- 

5t  rn>0  at  dt  dt  dt 

Equation  A.l  then  becomes 

dfl 

- =  e(Ej  +  v  x  Bx ) . v  f Q  (A. 3) 

dt 

We  now  have  only  one  variable,  t,  whereas  in  Euler's  coor¬ 
dinate  system,  we  would  have  three  independent  variables. 

The  perturbed  electron  distribution  function  is  then 


given  by: 
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f1{t)  =  e  |  [  E  -l  ( t '  )  +  v(t')  x  Bj  ( t '  )  J  .  vp  i  f  D  dt '  (A. 4) 

—  oo 


Along  the  unperturbed  tragectory,  r  Is  a  constant  and  the 
only  force  here  Is  the  magnetic  force  due  to  the  uniform  and 
constant  background  magnetic  field  BQzQ.  The  solution  is 
therefore  a  uniform  motion  along  z  and  a  circular  motion  in 
the  perpendicular  plane,  i.e. 


Vz  (t)  =  Vz  (t- ) 

Pt  (t)  =  pt  ( t ' )  (A. 5) 

r  (t)  =  r  (t1  ) 

Hence , 

3  a  3  3 

vp  =  (—  .  — )  =  v  ( —  '  — )  <A-6> 

3pt  <3pz  *  9pt’  apz’ 

and  therefore 


3fo 

3fo 

a  f  o 

- - 

xo  + 

— 

Yo  + 

—  z0 

3px 

apy 

3PZ 

3fo 

3pt 

3fo 

apt 

3fo 

- - 

— 

xo  + 

— 

Yo  + 

- z 

3pt 

<)px 

3P  t 

aPx 

apz 

3fo 

cos# 

3fo 

*fo 

= - 

xn  +  — 

sina  yQ 

+  — 

3pt 

3pt 

9pz 

(A. 7) 
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A. 2-1  Dispersion  Relation  for  TE  modes 

The  electromagnetic  wave  fields  can  be  described  by: 


E1  ■  Er  *0  +  E*  *o 

B1  -  Br  ro  +  B*  +  Bz  zo 


Then : 


E1  +  V  x  Bj  =  <Er  +  v0Bz  -  vzB0)  rQ  + 

<E0  +  vzBr  -  vrBz>  «o  +  <vrB0  "  v0Br>  zo  <A-8> 


For  a  TE  wave. 


B 


0 


B 


r 


(A. 9) 


Also , 


Vr  =  v^  COS  (0  -  (f)  ) 

=  Vt  sin (  0  -  0) 


(A. 10) 


Substituting  equations  A. 9  and  A. 10  into  equation  A. 8  we  get 


kv„ 


!1  +  V  x  Bj  =  [(1  - - )Er  +  vtBz  sin(0  -  0)j  rQ  + 


kv. 


(1  _  - )E^  _  vtBz  cos(0  -  0)]  0Q  + 
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Vtk  r 

- |Er  cos  (  9  -  0)  +  E  si  n(0 

U) 


(A. 11 ) 


Taking  the  dot  product  of  equation  A. 11  with  A. 7,  we  get 


kv_  ofn  vt<  c3fn 

(Ej  +  v  x  Ba)  .v  fQ  =  [(1 - ) - + - ] 

(0  3p^-  (D 

j^Er  cos(0  -  0)  +  sin(0  -  0)j  (A. 12) 

where  we  have  used  the  following  relationships: 

ro  •  xo  =  •  Yo  =  cos  * 

ro  ■  Vo  =  -  *  xo  =  sin  « 

(refer  to  figure  A.l  for  notation). 

Now,  the  solution  of  the  wave  equation  for  the  TE 
electric  field  is  given  by: 


m 


Er  = - Bzo  — -  Jm(ktr)  e 


j(wt  -  kz  -  m0) 


ktr 


t1 


(A. 13) 


J  ^zo  '  ^tr)  e 


j  ( to  t  -  kz  -  m0  ) 


With  the  help  of  Bessel’s  recursion  equations 
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ktr 


1 

2 


[Jm-l<ktr>  +'  Jm+l<ktr)] 


_  Jm+l(ktr)] 


and  Bessel's  addition  theorem 


-jm(0o  -  0) 

Jm^ktr^  =  e  -  Jm+p-ktro^  ^p^kt^L^ 

p=-°o 

K 

eJp<2  -  0O  +  0) 


where  we  have  used  the  geometrical  relationships 


,2  - 


2  2 
ro  +  rL 


2  rQ  rL  cos  a 


and  a  =  ^  -  0Q  +  0 

we  find  that 

Er  cos  (  9  -  0)  +  E0  sin  ( 6  -  0 )  =  j  ~  Bzo  e^wt  kZ '  e  J*m<^ 

x 

00  JP  ( “  "  </>0  + 

I  Jm+p(ktro)  Jp’(ktRL)  e  2  (A. 14) 

P  —  —  o° 


and  hence,  by  substituting  A.  14  into  A.  12  and  the  result 
into  equation  A. 3,  we  get: 
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dfl  ®B zo  9f0  3f0  j(»t  -  kz)  -jm0c 

- =  j - (w  -  kvz ) - +  vtk - 1  e  e 

dt  kt  3pt  3pz 


00  jp(r  -00  +  0) 

Z  J(ktr0)  Jp' (ktRL)  e  2 

P=-0O 


(A. 15) 


Note  that  the  only  time-dependent  functions  in  the  above 
expression  (equation  A. 15)  are 

71 

j(wt-kz)  jp(-  -  0O  -  0) 

e  and  e  2 

Integrating  equation  A. 15  with  respect  to  time  (from  t'=  -® 
to  t ) ,  we  get : 


fi(t)  = - e 


eBz0  j(wt  -  kz)  -jm<t>0  dfQ  3f0 

-  e  e  I ( m  -  kvz)  -  +  vtk  - 


Jm+p(ktro)  Jp ' ( ktRL  ^  -  «o  + 


Z  - - - - - e  2 

p=  — oo  C 0  -  kvz  +  p0o 


(A. 16) 


Equation  A.  16  represents  the  response  of  the  beam  to  the 
perturbation  introduced  by  the  wave  field. 

This  perturbation  in  the  electron  distribution  function 
gives  rise  to  a  perturbed  beam  current  which  will  be  the 
source  of  the  perturbation  in  the  wave  field.  Hence,  in  the 
wave  equation,  this  source  term  must  be  included  to  find  the 
response  of  the  field  to  the  perturbation  in  the  beam,  i.e.. 
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T 


(i 


nr 


Ym  2 

Kt  Rw 
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and  Jm,2(k^Rw)  =  0  because  of  the  boundary  conditions. 

Substituting  into  equation  A. 18,  we  get  the  following 


(k2  +  kt2) 


-j  (wt 
e 


where 


nr 


Kmn  -  <1  “ 


k  2  R  2 
Kt  Kw 


)V<ktRw) 


6J 

60 

(A. 19) 


The  integration  on  the  right  hand  side  needs  more  work, 
For  self  consistency,  the  current  perturbation  must  have  the 
same  time  dependance  as  the  field  that  created  it,  i.e. 

j(o>t  -  kz  -  m0) 

J  ^  =  J  ^  e 
r  ,  0  r  ,  0 

hence , 

aJr 

-  =  -jmJr 

60 


Also,  we  get  through  integration  by  parts: 
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/  “  t  )  Jm^  ktr ) dr  “  /  amnrJ0  Jm  '  ^  kt r 

«3r  r 

Hence,  the  terms  in  equation  A.  18  to  be  integrated  with 
respect  to  r  can  be  written  as: 

I  ["ktrJ0Jm'  (ktr>  +  jmJrJm<ktr>]dr  (A. 20) 

With  the  help  of  the  following  identities, 

Jm'(ktr)  =  -[Jjn-i  (ktr)  -  Jm+1(ktr)] 
m  1 

— Jm(ktr)  =  -[jmn(ktr)  +  Jm+1(ktr)] 
ktr  2 

the  integral  A. 20  becomes: 

-  j  Jr )  1  ^  amnr  )  ”  ^  ^  +  ^  Jr  ^m+1  ^  amnr  )  ]  ^  *  2 1 ) 

2  o 

Using 

Jr  “  -  efvrftd3p  and  J 0  =  -  efv^f^p 

and  also 

vr  =  vt  cos( 0  -  0)  and  sin(0  -  0) 


we  get : 
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Pt  j  (8  -  <!>) 

-  jJr  "  jej— e  flPtdptd0dpz 


Pt  -3(0  -  4>) 

Jd>  +  JJr  =  '  Jej  —  e  f tptdptd0dpz 

rm0 


Using  Bessel's  addition  theorem  : 


jm(0Q  -  0)  oo 

Jm(^tr^  =  €  £  *^m+p  r  l^tro 1  ^ t^L ^ 

P  '  =-CO 

K 

-jP! (-  -  0O  +  9) 
e  2 

we  find  that 

j(0~0)  - j ( 0  -  0 ) 

jm_i(ktr)e  +  Jm+i(ktr>e 


(A. 22) 


jm(0o  -  0)  » 


l  Jm+p'<kt ro! 


3  ’  p  *  (  )  e 


-  JP 1 (-  -  0O  +  0) 


(A. 23) 


Substituting  equations  A. 22  and  A. 23  into  A. 21,  and  the 
result  into  equation  A. 19,  we  come  up  with  : 


kD2  -  (k2  +  kt2)  *= 


,uekt  -j(ojt  -  kz  -  m0) 
2 - e 

Rw2Kmn 


rR w  r  r  r2n  pt  jm^o 

J  rdr  I  dpz  j  dpt  j  d0 - fje 

o  -«  o  o 


0) 
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co  -  jp '  ( -  -  0O  +  e ) 

£  Jm+P' (ktr)J'p- <ktRL)e  2 

*  ^  —  CO 


(A. 24) 


The  integration  with  respect  to  9  is  simple  since 
everything  is  constant  with  respect  to  9  except  for  the 
terms 


^  eJP(|  -  0O  + 

P  -oo 

which  are  in  the  expression  for  fj.  Hence,  we  have 


2/r  oo  =o  jp(_  -  0Q  +  9)  -jp'(-  -  <t>0  +  9) 

I  l  I  e  2  e  2 

o  P=~"  p'=-co 


j(p  -  P' ) (- 
=  l  e  2 

P.P’ 


0 0)  2,t  j(p  -  p’  )0 

J  e  d  8 


2k  j(p  -  p' )0 
Since  I  e 


delta , i . e . , 


PP' 


-  I? 


dO  =  2^(5ppt  where  <5ppi  is  the  kronecker 


P  t  P' 
P  =  P' 


we  get  upon  substitution  into  equation  (A. 24; 


oo  oo 


4 R.f  *>  oo  p 

kD2  -  (k2  +  kt2)  =  -  -----  J  rdr  /  dpz  /  dpt  — 


Rw^mn  o 


O  f  mo 


r,  ,afo  ,  «o,  Jm+p2(KtR0)Jp'2(ktRL) 

[  (a)  -  kv2) - +  vtk - J  Z - - - 

-)pt  p  u>  -  kv2  +  pfi 


(A. 25) 
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Let's  choose  the  unperturbed  distribution  function  to  be 


^  {  rO  ^O^^Pt'^Z^ 

2*r 


(A. 26) 


where  NQ  is  the  number  of  electrons  per  unit  length  of  the 
wave  guide,  the  Dirac  function  <5  ( rQ  -  RQ)  means  that  all  the 
electrons  have  the  same  guiding  center  located  on  a  cylind¬ 
rical  surface  defined  by  r  =  R0#  and  g(pt,pz)  is  an  arbit¬ 
rary  function  of  pz  satisfying  Jgd3  =  1. 

Then  (A. 25)  becomes  : 


kn2  -  (k2  +  kt2)  = 


2f/e^N0 


Pt 


^  Kmn  - 


j  dPz  /  dPt  —  [(•  "  kvz>-“-  +  vtkz - J 


r  mr 


3Pi 


-  Jm+p2(ktRo)Jp'2(ktRL) 

g(Pt<pz)  I  - 

p=~co  a)  -  kvz  +  pO 


ap2 

;  A . 27 ) 


Now,  we  need  to  integrate  the  right  hand  side,  i.e. 


“  Pt  r  kPZ  3  Pt  3  i 

J  dPz  J  dpt - [(“  - - ) - + - k - J 


rmQ  3pt  rmc 


ap. 


g(Pt<Pz)  l 

p=- 


Jm+p  2(ktR0)Jp' (ktRL) 


r® - P2  +  P°0 

nu 


(A. 28; 


where  we  have  let 


u)  -  kvz  +  pO 


r<D 


k 

—  Pz  +  P°o 

nu 


r 
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vz  =  pz/rm0 


and 


Vt  =  pt/rmc 


Integration  by  parts  of  (A. 28)  yields 


2jue^N0 


Pt 


kD2  -  (k2  +  kt2)  = - -  l  f  dpz  f  dpt  — 


V*mn  P~ 


00  —  oo 


r  mf 


[>■ 


2  (<02  -  k2c2  )  Hmp  ( x ,  y ) 


^vz  ^ Qmp^ x » Y ) 


( a  -  kv_  +  pO ) 


m  -  kv_  +  pfl 


]g(Pt  >Pj 


where 


y  =  ktRL 

Hmp<x-Y>  =  [Jp-m(x>V (Y)]2 

Qmp(x-y)  =  2HmP(x'Y)  +  y[Jp-m2(x)Jp'(y)Jp"(y)  + 

\  Jp-m-l2(x)Jp,(y)Jp-l"(Y)  - 

\  JP-m+i2(x)V (y)JP+i"(y)] 

Let  us  now  choose 

1 

g(pt,pz)  = - <5(pt,pto)<5(pz  -  pZQ) 

2*pt 


that  is,  the  momentum  distribution  function  of  the  electrons 
is  the  one  for  an  idealized  cold  electron  beam  (the  energy 
spread  is  very  small).  The  integration  is  then  easily  car¬ 
ried  out.  The  dispersion  relation  of  the  TE  modes  for  a 


gyrotron  is  then  : 
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B1  ~  Br  ro  +  B0 


Then 


v  7  <u 


E1  +  v  x  Bj  =  (1 - )(Er  rQ  +  E  4>0)  + 

c  kc 


and 


{e.,  +  —  —  [e^co s(0  -  0)  +  E  sin ( 0  -  0) 
1  '*  c  kc  1  r  * 


<E1  f  V  '<  Bl>  -vpfo  = 


vz  w  r  .  i  WJ“o 

(1  -  --  —  )|Ercos(0  -  0)  +  E  sin ( 0  -  0)j - + 

c  kc  ^  8p^ 


df. 


v  t  ^  u  L  o 

(ez  +  —  —  j^Ercos(0  -  0)  +  E  sin(0  -  0)jj - 

c  kc  ^  8pz 


But  now 


k  j(a>t  -  kz  -  m0) 

Er  *  -  EzoJm’<ktr>e 


m  k 

^  5=5  E2oJm^tr'e 


j(wt  -  kz  -  m0) 


kt  ktr 


After  substitution,  and  using  Bessel’s  addition  theorem 
we  get 


« 1 

dt 


eE, 


zo 


P  =  -o 


z  {  (k 


Rt  k 


Lt 


«'  vz  afo  vt  w'  afoi  , 

- )  —  +  ]  +  l} 

c  c  8pt  c  c  3p2 


j  ( «*t 

Jp(  ktRL.)  e 


K 


kz)  -jn»0o  jp(- 

J  (ktrQ)e  e  2 

m+p 


</>o  +  °) 


and,  upon  integration. 
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u)  -  kvz  +  pu 


(  A  .  3 1  ) 
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We  can  substitute 


N  „ue2 


Rw  Rmn 


by 


2c2K 


mn 


Integration  by  parts  of  equation  A. 31  yields 


where 


«  2 

4  v  k 

(-)  -  (k2 

+  kt2)  = - -  (-  j  -) 

c 

rKmnRw 

ffl  2 

[(j2  -  (k2  +  kt2)c2]GmN(x,y) 

[Pt 

(w  -  kvz  -  NO)2 

vz 

[w  - 

(k2  4-  kt2)— jRmN(x,y) 
k 

u>  -  k v 2  -  NQ 

X 

=  ktRo 

dzo 


B 


zo 


(A. 32) 


y  «  ktR 


t*L 


GmN ( x ' y ^  ”  jn(Y>jn'0) 


R, 


mN(x,y)  =  2GmN(x,y)  -  Jm_N2(x)  y  —  [-  JN(y)JN’(y) 


3 

ay 


A. 2-2  General  linear  dispersion  relation 

Combining  the  results  of  the  interaction  with  TE  modes 
(equation  A. 30)  and  the  result  of  the  interaction  with  TM 


microcopy  resolution  test  chart 

nation  At  BURCAU  or  sianoam*->9m-» 
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modes  (equation  A. 32),  we  finally  get  the  general  dispersion 
relation  : 


2  2 
-  -  ik*  *  kt*i  =  -  — - — - 
c  rKiTmc2 


{». 


jk  Ez°l 

r  2 

0) 

10  Bzo 

*vz  - 

NO) 

k  'in 

Z  O  r  r 

{<0  -  kvz)QNm(x,y)  -  j- - (0  -  (k2 

u>  B „„ 
zo 


a)  -  kv2  -  NO 


+  1c  2  )  -.?  j  R 


Nm 


} (A. 33) 


If  we  substitute  Ezo/Bzo  by  jcC/A,  we  then  get  a  real 
dispersion  relation 

Ezq  C  C* 

- =  -  *  j  C  - 

Bz0  A  A 

The  first  term  in  the  curled  brackets  on  the  RHS  repre¬ 
sents  the  contribution  from  the  instability,  the  second 
term  is  the  attenuation  due  to  the  Weibel  instability. 

A . 3  Simplification  of  the  dispersion  relation 

The  linear  theory  allows  us  to  derive  the  small  signal 
dispersion  relation  (eq.  A. 33)  for  all  types  of  beam-TE  and 
beam-TM  interactions.  It  may,  however  be  possible  to  make 
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some  approximations  in  order  to  have  a  simpler  equation  to 
work  with.  To  achieve  this  purpose,  we  should  compare  the 
contributions  from  each  interaction  to  the  one  which  is 
dominant  (i.e.,  the  TE^-beam  interaction).  We  know  that  the 
magnitudes  Ap/AQ  and  Cp 1 / kQ  are  small  numbers  since  only  the 
fundamental  mode  is  propagating,  the  higher  order  modes 
being  small  and  evanescent. 

However,  this  is  not  enough  to  neglect  the  TM  terms. 
Indeed,  we  need  to  know  what  is  the  contribution  of  the  TM 
modes  to  the  coupling.  To  do  this,  let's  calculate  the 
ratios  of  the  TM  to  the  TE  contributions  for  each  term  on 
the  right  hand  side  of  the  dispersion  relation  (equation 
A. 33)  : 

TMI  2k2  Ak  GNm  C  ' 

TE I  Gain  0t2  kQ  kQ(l  -  k2/kQ2)  HNm  AQ  . 

TM!  k  1  -  0Z  k0/k  -  20zAk  RNm  C  ' 

Ew  =  * - ? 

TElWeibel  kQ  (1  -  k/ko)0z  QNm  AQ 

For  a  given  (^,k),  these  contributions  are  negligibly 
small.  We  can  therefore  disregard  the  contribution  of  the  TM 
Modes  to  the  dispersion  relation,  the  error  thus  introduced 
being  small. 
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APPENDIX  B 

DERIVATION  OF  THE  NON-LINEAR 
EVOLUTION  OF  THE  RADIATION  FIELDS 

B ■ 1  General  background : 

The  Investigation  of  the  non-linear  evolution  of  the 
CMI  starts  with  solving  the  equation  of  the  single  electron 
motion  in  the  radiation  fields.  It  is  assumed  that  the  phase 
trapping  of  the  bunching  electrons  in  the  wave  fields  is  the 
mechanism  responsible  for  the  saturation.  The  relations  that 
determine  the  collective  response  of  the  electrons  to  the 
radiation  fields  are  derived  following  the  averaging  of  the 
results  over  the  initial  random  phase  distribution  of  the 
electron  transverse  velocity.  The  response  of  the  radiation 
fields  to  the  current  induced  by  the  electron  bunching  is 
then  analyzed  using  the  equation  of  conservation  of  energy. 
This  self-consistant  approach  leads  to  the  derivation  of  a 
single  non-linear  equation  that  describes  the  evolution  of 
the  radiation  field  intensity. 

B.2  Characteristic  equations: 

The  motion  of  the  electrons  in  an  electromagnetic  field 
superimposed  on  a  d.c.  magnetic  field  is  governed  by  the 
following  set  of  coupled  non-linear  equations 
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where  is  the  n-th  root  of  the  equation  Jm,(k<tRw)  =  0. 

Referring  to  figure  A, 1  for  notation,  the  spatial  coor¬ 
dinates  (R ,<t>)  of  the  electron  can  be  expressed  in  terms  of 
the  coordinates  (RO,0O)  of  the  guiding  center  of  the  elec- 


93 


iron  and  the  coordinates  (R^,#)  of  the  motion  of  the  elec¬ 
tron  about  the  guiding  center,  i.e., 

r2  =  K02  +  rL2  ~  2RoRLcos<2  '  +  ») 

0  =  0O  -  Sin-1  (Rl/Ro)cos(0o  -  0) 

whore  R^  =  vt/0  =  pt/m0r0  is  the  Larmor  radius  of  the  elec¬ 
tron  . 

If  the  wave  frequency  m  is  near  one  of  the  Doppler 
shifted  harmonics,  i.e.,  w  =  kvz  +  Nfi ,  the  slowest  term  of 
the  forces  experienced  by  the  gyrating  electrons  has  the 
dominant  influence  on  the  orbit  evolution  of  the  electrons. 
Retaining  this  term  only  in  equation  B.l,  we  first  define  a 
self  consistent  trajectory 

rl  =  rl^) 

t 

0  —  0  ( t )  =  0 0  +  0(t)  +|  Q ( t 1  ) dt 1 

o 

prior  to  express  B.l  explicitely. 

RL<t)  and  0(t)  are  to  be  determined  self-consistent ly 
and  represent  respectively  the  time  dependant  Larmour  radius 
and  that  term  in  the  phase  that  varies  slowly  in  response  to 
the  slow  time  varying  force. 

We  then  substitute  these  relations  and  the  expressions 
for  the  fields  ( eq .  B.4)  into  B.l  and  B.2  and  neglect  all 

the  fast  oscillating  terms. 
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d  0o2a»  k2c2 

--vz  =  -  a - ( kvz  - )«JN'(a)cos0  (B.9) 

dt  r2kkt2c2  ui 


where 


a  =  kt.RL 


eBo 

°o  =  TO  =  — 


m0c 


t  UK 

0  =  (m  -  N)0O  -  [  Awft'Jdt1  +  N (0  +  0Q  +  -)  +  kzQ  +  - 
o  2  2 


r  -  rf 


Aa>  =  m  -  NO  -  kvz  =  Aa>0  +  ( a)Q  - 


Aw0  =  ce  -  wQ 


N0o 

"o  =  f  kVZO 


a  ""  (  ~  ^m-N^t^o- 


B.9 


Under  the  grazing  condition,  v  =  ko^/w  =  vn ,  equation 


dv. 


reduces  to  -  =  0  and  thus  vy  -  vyo 

d  t 


To  close  the  system  of  equations,  it  is  necessary  to 
include  a  field  equation  that  takes  into  account  the  dynami  - 
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cal  effects  of  the  medium  on  td\e  fields.  The  energy  conser¬ 
vation  equation  B.3  is  used  for  this  purpose,  and  under 
grazing  conditions  it  becomes 

dr  o e o  d 

n0m0c2—  + - -  Eq2  =  0  (B.10) 

dt  2  dt 


nQ  is  the  electron  beam  density  averaged  over  the 
cross-section  of  the  guide. 

<  >  stands  for  an  average  over  the  random  initial  phase 
distribution  of  the  electron  beam. 
w 

Eo  ”  Bzo 

k^c 


m4, 

«  *  (' - ;)J,„2(KtRw»  Is 

“t  Rw 

represents  the  energy  distribution  ove 
the  guide. 

The  integration  of  equation  B.10 


filling  factor  and 

r  the  crr'sc  -section 


leads  to 


of 


oe. 


n0m0c‘ 


( <r>  -  < 


O  r 

rQ>)  = - 11 


( t) 


(0) 


B  .  1 1 


B.3  derivation  of  the  NL  field  equation: 

From  the  equation  of  conservation  of  energy,  we  can 
find  the  polarization  current 


density ,  i . e . , 
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(B. 14) 
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We  assume  in  the  present  analysis  that  the  effect  of 
depletion  of  the  electron  rotational  free  energy  to  be 
insignif icant .  Then,  EQ ( t )  and  0Q( t)  are  expected  to  vary 
with  time  much  faster  than  the  other  variables  a,  y,  and  v2 . 
We  can  then  neglect  the  slow  time  variation  of  a  and  y  and 
make  a  *  aQ  in  equation  B.6  and  7  =:  7Q  in  equation  B.8. 

The  averaged  (ever  the  random  initial  phase  distribu¬ 
tion)  polarization  current  density  will  then  be 


<jp> 


-  DN - ~-Jm-N<ktRo  )  <- Jjj  r  (  a  )  cos (/>> 


en^Q 


00 


<jp>  =  - 


i) 


N 


enoQo 


'm-N  (  ktRo  )  _~JN  '  («0  )  <coS(/>>  {  E  .  15) 


Using  equation  B.14  and  writing  equation  B.15  as 


<Jp>  -  -  (-1)NP<COS0> 


we  get 


<COS0> 


<7£o  dEo 
(-  1 )NP  dt 


<COS0>  =  (-1) 


N+l 


(T6C  dE0 


P  dt 


:b. i6) 


where 


enono  «o 

,*  Jm-N<ktRo^“  'Jn'  (ao' 
kt  rQ 
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Equation  B.8,  with  the  simplifying  assumptions  a  =  aQ 

and 

r  -  rQ.  can  now  be  integrated  to  obtain 


a 


r  -  r0  -  -\'(«0)(  dt  'a(t  ’  }cos0  1  (H.17) 


k  2_  2 
Kt  c  ro 


We  can  also  approximate  the  expression  for  Ao»  (using  r 
-  70)  to  give 

7  "  ^o 


k2c2 


Aa>  =  AwQ  + 


>o - ) 

w 


( B  .  1 8  ) 


which,  upon  substituting  equation  B.17,  yields 


A  u)  ^  A  a) 


0  +  A QJ  dt'a'cos0 


where 


Aq  =  2  (u> 2  -  k2c2)a0JN'  («D) 

roktc 


a'  ■  a(t ’ ) 


0'  «  0(f) 

Using  equations  B.16  and  B,17,  the  average  of  Aw  can  be 
expressed  as 


<7£, 


<Aw>  =  A - 


w  2  J'onomoc' 


[ED2(t)  -  E02 (0) 

(b. iq: 


100 


Taking  into  account  only  the  resonant  processes,  the 
averaged  polarization  current  density  will  be  shown  to  obey 
a  second  order  differential  equation. 


From  the  definition  of  Jp,  we  get 


d2  a 


d2jp  f 

----  =  p 

rd  <t> 

dt2  1 

dt 

<«)] 

’t2  r 

«/7  JN'(a) 


j  cos <p 


[ 


+  2 


‘dt 


d  a 

"[~jn' (<*)] 

S<fiv  dt  r 

dt  a/r  (  a  ) 


j  sin0  J 


d2  a 


d2 


dt  * 


d  02 

dt2*; 

r 

'  («)] 

P  \  < - COS0>  - 

< - 

— 

1  dt 

a/r 

JN 

'  («) 

d  a 

d20  d0 

dt 

jn' 

(«)] 

- sin0>  +-  2 <  — 

— 

— 

- s 

dt 2  dt 

a/r 

JN 

1  (a) 

cos</>>  + 


in0>j  (B. 


20; 


In  the  following,  we  are  going  to  derive  the  explicit 
expressions  of  functions  averaged  over  the  random  initial 
phase  distribution  of  the  electron  beam.  Only  the  linear 
terms,  and  the  non-linear  terms  of  lowest  order  are  retai¬ 
ned  . 

Using  equation  B.6  through  B.8,  we  get 
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We  also  have 
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d0 

—  =  -  A  a) 
dt 


[l  +  aJN(a)sintf>]  -  aDQsin0  (B.23) 


-  =  (Aw0)^  +  Aq ( Am  +  Aw0 ) 

dt 


j  dt '  a  '  cos<t> 1  +  a2DQ2sin20  +  2aDoAaisin0  (B.24) 


d20  da 

—  -  =  -  aAQcos 0  -  DQ ( — sin0  -  aAwcos </>)  + 


dt' 


dt 


— [Do2  -  Dl2  "  AoJN(«)]si 


sin2<x 


(B.25) 


Then, 


d0  ^ 

< - cos 0 >  =  (  A0io)  2<COS0>  +  a2DQ2  <sin20cos0>  + 

dt 

aD  <Ao)sin20>  +  An<(Aw  +  Aw  )cos0 


j  dt'a'cos 0’> 


(  B  .  2  6  ) 


d20  1  1  da  1 

< - sin0>  =  -  -aAo<sin20>  -  -  — DQ  +  -aDQ<Awsin20> 

dt2  2  2  dt  2 


+  -a2(D02  -  Dj2)<cos0> 

4 


(B. 27) 


The  following  simplifications 


<sin20> 


0 
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2  1  1  1 

<sln^0cos0>  -  -<cos 0>  -  -<cos30>  =  -<cos0> 

4  4  4 

3  3  1  3 

<cosJ0>  =  -<COS0>  +  -<cos30>  =  -<cos0> 

4  4  4 


can  be  substituted  into  equations  B.21,  B.22,  B.26  and  B.27, 
and  the  result  into  equation  B.20  which  then  becomes: 


d  JP  r  r  2  Aow  3  Aow 

dt2  1  L  *2  _  k2c2  o  4  _  k2c2 


0  7  0  0*1  ^  AoW 

‘  4.  1  /A„  \Z  _ _ /  n 


(Dq2  +  D^Ja2  +  (Ao>0)2] - (Dc 


2  (o2  -  k2c2 


3  Aow 

-<Aojsin20>a - - + 

4  (o2  -  k2c2 


V. 

Aq<  { Aw  4-  Aa>o)cos0j  dtfafcos0'>  (B.28) 


We  still  need  to  find  the  explicit  expressions  for 


*  7  c L 

t>sin  0  >  and  <(Aw  +  Ag>q)cos0  J  dt,a,cos0'> 


k. 

<Aa»sin20>  =  AQ< 2sin0cos0|  dt'a'cos0'> 


L 

-  AQj<sin0j  dt ' a ' cos0cos0 1 >  + 
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These  approximations  are  legitimate  if  only  the  linear  and 
the  lowest  order  non-linear  terms  are  taken  into  considera¬ 
tion. 

With  the  help  of  these  approximations,  we  get 

1  t 

<Awsin20>  -  -AQ[ <sin0> |  dt'aTcos<0  -  0f>  + 

2  1  o 

t 

<cos</)>J  dt,a,sin<0  -  0’>J 
o 
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1  t 

=  -  (  <  Aa>>  +  Ao)Q)  f  dt,a,cos<0  -  0  1  > 

2  o 


Also,  from  the  definition  of  Jp,  we  have 


dJ  enQ0o 

—  - - Jm-N<ktRo> 

dt  kt 


d  a  a  d0 

| cos0 — -  -JjfMa) — sin^l 

dt L  r  r  dt 


dJp  enoflo 

<  >  _  _Jm~N^ktRo^ 

dt  kt 


a  d  0  d  a 

[<-JM'(a) — sin0>  -  <cos0 — [-JN,(a)>[J 
r  dt  dt  r 


dJ_  d0  1  i 

<  —  ->  =  Pj  < — sin0>  -  -a(D0  -  - - ---)  ( B  .  29  ) 

dt  L  dt  2  w2  -  kzc2 


Since 


d0 

^  op 

--sin0  =  -  Awsin0  -  aAwJN( a ) sin^0  -  aDQsin^0 
dt 


We  have 


d  10 

<--sin0>  =  -  <A<osin0>  -  -aDQ 
dt  2 

1  t  1 

=  -  Aa>o<sin0>  -  -Aq[  dt,a,sin<0  -  0’>  -  -aDc 

2  o  2 
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which  yields 


d  0 

< — sin  0> 


dt  1  Aq  t  1  aD0 

<sin0>  =  - - -  - - j  dt  '  a '  s  in<  0  -  0'>  -  - - 

Aa>0  2  AwQ  o  2Aw0 


(B. 30) 

Using  equations  B.29  and  B.30,  we  get  the  following 
expression 


<sin0> 


A0(« 

<u2  -  k2c2 


1 

2  A“o 


1  Ao 

-  -  J  dt'a!sin<0  -  0'>  (B.31) 

2  Aoiq  o 


We  also  have  from  equation  B.13 


<COS0> 


(  B  .  3  2  ) 


Substituting  equations  B.31  and  B.32  into  equation 
B . 28 ,  we  get 


d2  J  3  a  A02a>  dJ 

< - ->  +  - - 1_<  —  -> 

dt2  4  A<oq  at2  -  k2c2  dt 


3  A0‘(« 

+  -a - 1*. 

4  <o2  -  k2c2 
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AQw  3  A  quj 

-  a2(DQ2  +  DX2  )  +  3a2  - - -—(D,,  ~  "  "i - r5”5>]JP 

at*  -  k*c*  4  ai*  -  k*c* 

3  A  O2o>  P  Aqw 

-a~-  ~ -  ~  1(^0  ”  “  ~  ”  ~  0 )  +  A0IS  +  a^oj^c  + 

8  u)2  -  kzcz  ^u)0  u)Z  -  kzcz 

1  da  1  AQw 

-PA0(ao>0  +  <Ao»)ic  -  P~(D0 - r - ~r)  (B.33) 

2  dt  2  wz  -  kzc^ 


where 


t 

Ic  =  |  dt,a,cos<0  -  1  > 

o 

t 

Is  =  |  dt'a'sin <0  -  0’> 
o 


If  we  let 


C 


o 


> 


P(D0 


1  Aqw 

2  <u2  -  k2c2 


2S 


3  aA02a» 

4  A<u0(w2  -  k2c2) 


Aw 


1*  =  Ao,02  -  a2  [d j 2  +  (Dc 


A0io 


2d)2  -  k2c2 


>2] 


3  Aq2(u 

-a - ---I,, 

4  <a2  -  k2c2 
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equation  B.33  becomes 


dJ, 


p  ~~p  o 

< - £>  +  2S  —  —  +  (  A<o1  )  ZJ 

dt2  dt  y 


da 


A0w 


-  CD-  -  SP[aDc  +  (Dq  -  + 

dt  m2  -  k2c2 


AoXs]  +  -PA0Ic(<A.>  +  Aw0] 


;B.34) 


Equation  B.14  relates  the  polarization  current  density 
and  the  radiation  field  intensity.  Using  this,  and  the  fact 
that  A  a)  plays  a  role  similar  to  the  population  inversion 
function  in  a  two-level  system  which  allows  us  to  define  a 
population  inversion  function 

nO  Y omoc 

W  = - (<Au>  +  A«i  )  (  B  .  35  ) 

2H<«02 

where  is  Planck's  constant. 

We  finally  get  a  single  equation  that  governs  the  non¬ 
linear  evolution  of  the  radiation  field  intensity 


(A<u12  +  C,  )  —  ]eo  = 

dt 


no 


p 


111 

id, 

D0  =  (w  -  kvzo)  — - |aoJN,^ftoM 

a°  da°  “>  2\ 

AO.J2  =  Aa-02  -  a2{D12  +  [dq  -  ^ 

3  w 

-  a  ( A02  "5 - 2_2)Is 

4  (u2  -  k^c*1 

a  =  (-  DN  -WktRo) 

Bo 

D,2  -  <•  -  kvzo)2  JK'(«0)---  I"  ---[«oJ»'<«ol 

1  daQ  aQ  daQ 

3  aAD2w 

S  *  "  ry  TO,  C 

8  Auiotw'1  -  k^c^ ) 

t 

Xc  =  J  dt 1  a  ( t f  )  cos< A0 ( t  -  t’)> 
o 

t 

Is  =  |  dt 1  a(t!)  sin<A0(t  -  t')> 

o 

m2  2 

a  =  (1  - 

kt2Rw2 


..{ . 


112 


BIBLIOGRAPHY 


[  1 ]  The  gyrotron 

V.A.  Flyagin,  A.V.  Graponov,  M.I.  Peletin,  and 
V.K.  Yulpatov. 

IEEE-MTT ,  Vol .  MTT-25,  #6,  June  1977 


[2]  The  Electron  cyclotron  maser.  An  historical  survey 
J.L,  Hirschfield  and  V.L.  Granatstein, 

IEEE-MTT,  Vol.  MTT-25,  #6,  June  1977 


[3]  Comparative  study  of  the  axial  and  azimuthal 
bunching  mechanisms  in  e 1 ec t romagne t  i c  cyclotron 
instabilities 

K.R.  Chu  and  J.L.  Hirschfield. 

Physics  of  Fluids,  Vol.  21,  #3,  March  1978 


[4]  Simple  macroscopic  theory  of  cyclotron  maser 
ins tabi 1 i t ies 
Y.Y.  Lau. 

IEEE-Elect ron  Devices,  Vol  ED-29,  #2,  February 

1982 


[5]  A  unified  theory  of  diocotron,  cyclotron  maser, 
and  negative-mass  instabilities 
Y.Y.  Lau. 

lEEE-Electron  Devices,  Vol.  ED-31  #3,  March  1984 


[6]  Relativistic  effect  in  electron  c yclotron  trans¬ 
verse  wave  devices 
G.J.  Sehn  and  R.E.  Hayes. 
lEEE-Electron  Devices,  December  1969 


[7]  Enhanced  microwave  emission  due  to  the  transverse 
energy  of  a  relativistic  electron  beam 
M.  Friedman,  D.A.  Hammer,  W.M.  Manhelmer  and 
P.S.  Sprangle. 

Physical  Review  Letters,  Vol.  31  #12,  September 

1973 


113 


[8]  Stimulated  emission  of  radiation* by  relativistic 
electrons  in  a  magnetic  field 
J.  Schneider. 

Physical  Review  Letters,  Vol.  2  #12,  June  1959 


[9]  Radiation  transfer  and  the  possibility  of  negative 
absorbtion  in  radio  astronomy 
R.Q.  Twiss. 

Austrian  Journal  of  Physics,  Vol.  11  pp .  564-579, 

1958 


[10]  The  gyrotron  :  high  power  source  of  millimeter  and 
submillimeter  waves 

A. A.  Andronov,  V.A.  Flyagin,  A.V.  Gaponov, 

A.L.  Goldenberg,  M.I.  Peletin,  V.G.  Osof,  and 
V.K  Yulpatov. 

Infrared  Physics,  Vol  18  pp .  385-395,  1978 


[11]  Cyclotron  autoresonance  masers  of  millimeter-wave 
range 

V.L.  Bratman,  G.G.  Denisov,  and  M.M.  Ofitserov. 
Relativistic  H.F.  Electronics  (Institute  of  Ap¬ 
plied  Physics,  Gorky),  1983 


[12]  Spatial  and  temporal  coherence  of  a  35  GHz  gyro- 
monotron  using  the  TE02  circular  mode 
M.E.  Read,  R.M.  Gilgenbach,  R.F.  Lucay  Jr.,  K.R. 
Chu,  A . T .  Drobot,  and  V.L.  Granatstein. 

IEEE-MTT ,  Vol.  28  #8,  August  1980 


[13]  A  high  power  gyrotron  operating  in  the  TEo41  mode 
B.  Arfin,  K.R.  Chu,-  D.  Dialetis,  and  M.'E.  Read. 
IEEE-Elec tron  Devices,  Vol.  ED-29  #12,  december 

1982 


[14]  A  Gyrotron  study  program 

G.  Boucher,  P.  Boulanger,  P.  Charbit,  G.  Fail Ion, 
A.  Herscovici,  E.  Kammerer,  and  G.  Mourier .  In¬ 
frared  and  Millimeter  waves,  Vol.  9,  1983 


114 


[15]  The  problems  of  increase  in  power,  efficiency  and 
frequency  of  gyrotrons  for  plasma  investigation 
A . S .  Fix,  V . A .  Flyagin,  A .  L .  Goldenberg, 

V.I.  Khizhnyak,  S.A.  Malygin,  S.E.  Tsimring,  and 
V.E.  Zapevalov. 

International  Journal  of  Electronics,  Vol.  57  #6, 
1984 


[16]  Gyrotrons  for  high  power  millimeter  wave  genera¬ 
tion 

H.R.  Jory,  F.  Friedlander,  S.J.  Hegji, 

J.F.  Shively,  and  R.S.  Symons. 


[17]  Cyclotron  resonance  devices 
R.S.  Symons,  and  H.R.  Jory. 

Advances  in  Electronics  and  ElectroPhysics ,  Vol. 
55,  1981 


[18]  Gyrotron-TWT  operating  characteristics 

P.E.  Ferguson,  G.  Valier,  and  R.S.  Symons. 
IEEE-MTT,  Vol  MTT-29  #8,  August  1981 


[19]  An  experimental  gyro-TWT 

R.S.  Symons,  H.R.  Jory,  S.J.  Hegji,  and  P.E.  Fergtison. 
IEEE-MTT,  Vol.  MTT-29  #3,  March  1981 


[20]  Measurements  of  gain  for  slow  cyclotron  waves  on  an 
annular  electron  beam 

H.  Guo,  L.  Chen,  H.  Keren,  J.L.  Hirschfield,  S.Y.  Park, 
and  K.R.  Chu. 

Physical  Review  Letters,  Vol.  49  #10,  September  1982 


[21]  An  experimental  wide-band  gyrotron  traveling  wave 
amplifier 

L.R.  Barnett,  Y.Y.  Lau ,  K.R.  Chu.  and  V.L.  Granatsteiri . 
IEEE-Electron  Devices,  Vol.  ED-23  #7,  July  1981 


[22]  Theory  of  a  wide-band  distributed  gyrotron  traveling 
wave  amplifier 

K.R.  Chu,  Y.Y.  Lau,  L.R.  Barrett,  and  V.L.  Granatstein. 
IEEE-Electron  Devices,  Vol.  ED-28  #7,  July  1981 


115 


[23]  Large  Signal  theory  of  a  two-stage  wide  band  gyro-TWT 
A.K.  Ganguly  and  S.  Ahn . 

IEEE-Electron  Devices,  Vol .  ED-31  #4,  April  1984 


[24]  An  active  circulator-gyrotron  traveling  wave  amplifier 
Y.  Y.  Lau,  L.R.  Barnett,  and  J.M.  Baird 
IEEE-Electron  Devices,  Vol.  ED-31  #3,  March  1984 


[25]  The  helically  wrapped  circular  waveguide 
K.J.  Bunch  and  R.W.  Grow. 

IEEE-Electron  Devices,  Vol.  ED-34  #8,  August  1987 


[26]  Theory  of  gyrotron  amplifiers  in  disc  or  helix  loaded 
waveguides 

J.Y  Choe  and  H.S.  Uhm . 

International  Journal  of  Electronics,  Vol.  53  #6,  1982 


[27]  The  theory  of  disc-loaded  waveguides 
E . L .  Chu  and  W.W.  Hansen. 

Journal  of  Applied  Physics,  Vol.  18,  November  1947 


[28]  Multimode  oscillation  and  mode  competition  in  high 
frequency  gyrotrons 

K.E.  Kreischer,  R.J.  Temkin,  H.R.  Fetterman,  and 
W.J.  Mulligan. 

IEEE-MTT ,  Vol.  MTT-32  #5,  May  1984 


[29]  Analysis  of  the  electron  cyclotron  maser  instability 
S.P.  Kuo  and  B.R.  Cheo . 

Physics  Letters,  Vol.  103A  #9,  July  1984 


[30]  Analysis  of  the  harmonic  gyrotron  traveling  wave  ampli¬ 
fier 

S.P.  Kuo,  S.C.  Kuo,  B.R.  Cheo,  and  M.C.  Lee. 
International  Journal  of  Infrared  and  Millimeters 
Waves,  Vol.  7  #4,  1986 


[31]  Theory  and  simulation  of  the  gyrotron  traveling  wave 
amplifier  operating  at  cyclotron  harmonics 
K.R.  Chu,  A.T.  Drobot,  H.H.  Szu,  and  P.  Sprangle. 
IEEE-MTT,  Vol.  MTT-28  #4,  April  1980 


116 

[32]  Waveguide  handbook 

Edited  by  N.  Marcuvitz. 

McGraw  Hill ,  1951 


[33]  The  reflection  of  an  electromagnetic  plane  wave  by  an 
infinite  set  of  plates 
J.F.  Carlson  and  A.E.  Heins. 

Quarterly  of  Applied  Mathematics,  Vol  4  #4,  January 

1947 


[34]  Self-consistent  large  signal  theory  of  the  gyrotron 
traveling  wave  amplifier 
A.K.  Ganguly  and  S.  Ahn. 

International  Journal  of  Electronics,  Vol.  53  #6,  1982 


[35]  Electromagnetic  slow  wave  systems 
P.M.  Bevensee. 

Wiley,  1964 


[36]  General  mode  analysis  of  a  gyrotron  dispersion  equation 
J.Y.  Choe  and  S.  Ahn. 

IEEE-Electron  Devices,  Vol.  ED-28  #1,  Jnuary  1981 


[37]  The  linear  and  self  consistent  non-linear  theory  of  the 
electron  cyclotron  maser  instability 
P.  Sprangle  and  A.T.  Drobot. 

IEEE-MTT ,  Vol.  MTT-25  #6,  June  1977 


[38]  Coherent  non-linear  theory  of  a  cyclotron  instability 
P.  Sprangle  and  W.  M.  Manheimer. 

Physics  of  fluids,  Vol.  18  #2,  February . 1 975 


[39]  A  note  on  gyrotron  traveling  wave  amplifiers  using 
rectangular  waveguides 
Y.Y.  Lau  and  L.R.  Barnett. 

IEEE-Electron  Devices,  Vol.  ED-30  #8,  August  1983 


[40]  AC  Space  charge  effects  in  gyrotrons 
G.  Dohler. 

International  Journal  of  Electronics,  Vol.  56  #4,  1984 


117 


[41}  The  lasing  mechanism  in  fast-wave  devices 
G.  Dohler. 

IEEE-Electron  Devices,  Vol .  ED-31  #3,  March  1984 


[42]  Measured  performance  of  gyrotron  oscillators  and  ampli¬ 
fiers 

V.L.  Granatstein,  M.E.  Read,  and  L.  Barnett. 

Infrared  and  Millimeter  Waves,  Vol.  5,  1982 


[43]  Status  review  of  research  on  millimeter  wave  tubes 
J.  Feinstein  and  K.  Felch. 

IEEE-Electron  Devices,  Vol.  ED-34  #2,  February  1987 


[44]  Study  of  noise  amplification  mechanism  in  Gyrotrons 
K.R.  Chen  and  K.R.  Chu. 

IEEE-MTT ,  Vol.  MTT-34  #1,  January  1986 


[45]  Field  theory  of  guided  waves 
R.  Collins. 

McGraw  Hill 


[46]  A  flush-mounted  leaky-wave  antenna  with  predictable 
patterns 
R . C .  Honey . 

I.R.E.  Trans,  on  Antennas  and  Propagation,  October  1959 


[47]  Theory  and  numerical  modeling  of  a  compact  low  field 
high  frequency  gyrotron 

P.  Vitello,  W.H.  Miner,  and  A . T .  Drobot. 

IEEE-MTT,  Vol.  MTT-32  #4,  April  1984 


[48]  Realization  of  a  stable  and  highly  efficient  gyrotron 
for  controlled  fusion  research 
Y.  Carmel,  K.R.  Chu,  M.  Read,  A.K.  Ganguly, 

D.  Dialetis,*R.  Seeley,  J.S.  Levine,  and 
V.L.  Granatstein. 

Physical  Review  Letters,  Vol.  50  #2,  January  1983 


[49]  Multimode  analysis  of  quasi-opt ical  gyrotrons  and 
gyroklystrons 

A.  Bondeson,  W.M.  Manheimer,  and  E.  Ott. 

Infrared  and  Millimeter  Waves,  Vol.  9 


118 


[50]  Practical  considerations  in  the  'design  of  a  higher 
power  1mm  gyromonotron 

J.D.  Silverstein,  M.E.  Read,  K.R.  Chu,  and  A.T.  Drobot . 
IEEE-MTT,  Vol .  MTT-28  #9,  September  1980 


[51]  A  125-260-GHz  gyrotron 

G .  F .  Brand-,  N.G.  Douglas,  M.  Gross,  J.Y.L.  Ma, 
L.C.  Robinson,  and  C.  Zhiyi. 

IEEE-MTT,  Vol.  MTT-32  #1,  January  1984 


[52]  Methods  of  efficiency  enhancement  and  scaling  for  the 
gyrotron  oscillator 

K.R.  Chu ,  M . E .  Read ,  and  A . K .  Gangu 1 y . 

IEEE-MTT,  Vol.  MTT-28  #4,  April  1980 


[53]  Experimental  examination  of  the  enhancement  of  gyrotron 
efficiencies  by  use  of  profiled  magnetic  fields 
M.E.  Read,  K.R.  Chu,  and  A.J.  Dudas . 

IEEE-MTT,  Vol.  MTT-30  #1,  January  1982 


[54]  Loaded  Q's  and  field  profiles  of  tapered  axisymetric 
gyrotron  cavities 

H.  Derfler,  T.J.  Grant,  and  D.S.  Stone. 

IEEE-Electron  Devices,  Vol.  ED-29  #12,  December  1982 


